o 
o 

(N 

< 






I 

o 



O 

o 

B 

o 
o 



X 



Angle-Resolved Photoemission Spectroscopy on Electronic Structure and 
Electron-Phonon Coupling in Cuprate Superconductors 

X. J. Zhou/'2,3 J c^i^ 1 J Devereaux,'^ N. Nagaosa,^ and Z.-X. Shen^ 

^ Dept. of Physics, Applied Physics and Stanford Synchrotron Radiation Laboratory, Stanford University, Stanford, 

CA 94305 

^Advanced Light Source, Lawrence Berkeley National Lab, Berkeley, CA 94720 

'^State Key Laboratory for Superconductivity, Beijing National Laboratory for Condensed Matter Physics and Institute of 

Physics, Chinese Academy of Sciences, Beijing 100080, China 

'' Department of Physics, University of Waterloo, Ontario, Canada N2L 3GI 

^ CREST, Department of Applied Physics, University of Tokyo, Bunkyo-ku, Tokyo 113-8656, Japan 






Contents 

I. Introduction 

II. Angle-Resolved Photoemission Spectroscopy 

A. Principle 

B. Technique 

III. Electronic Structures of High Temperature 
Superconductors 

A. Basic Crystal Structure and Electronic Structure 

B. Brief Summary of Some Latest ARPES Results 

IV. Electron-Phonon Coupling in High Temperature 
Superconductors 

A. Brief Survey of Electron-Phonon Coupling in 
High- Temperature Superconductors 

B. Electron-Phonon Coupling: Theory 

1 . General 

2. Weak Coupling - Perturbative and Self-Energy 
Description 

3. Strong Coupling - Polaron 

C. Band Renormalization and Quasiparticle Lifetime 
Effects 

1. El-Ph CoupUng Along the (0,0)-(7r,7r) Nodal 
Direction 

2. Multiple Modes in the Electron Self-Energy 

3. El-Ph Coupling Near the (7r,0) Antinodal Region 

4. Anisotropic El-Ph Coupling 

D. Polaronic Behavior 

1. Polaronic Behavior in Parent Compounds 

2. Doping Dependence: From Z~0 Polaron to Finite '. 
Quasiparticles 

3. Doping Evolution of Fermi Surface: 
Nodal-Antinodal Dichotomy 

E. Electron-Phonon Coupling and High Temperature 
Superconductivity 

V. Summary 

Acknowledgments 

References 



I. INTRODUCTION 

In addition to the record high superconducting tran- 
sition temperature (Tc), high temperature cuprate 
superconductors! l| |2) are characterized by their un- 
usual superconducting properties below Tc, and anoma- 
lous normal state properties above Tc. In the super- 



conducting state, although it has long been realized 
that superconductivity still involves Cooper pairs (^, 
as in the traditional BCS theory(|2 |a_^la), the ex- 
perimentally determined d— wave pairingrr) is different 
from the usual s— wave pairing found in conventional 
superconductors^ Mj The identification of the pair- 
ing mechanism in cuprate superconductors remains an 
outstanding issue IJIOJ ). The normal state properties, par- 
ticularly in the underdoped region, have been found to 
be at odd with conventional metals which is usually de- 
scribed by Fermi liquid theory; instead, the normal state 
at optimal doping fits better with the marginal Fermi liq- 
uid phenomenology yTj). Most notable is the observation 
of t he p seudogap state in the underdoped region above 
Tc (|13) . As in other strongly correlated electrons sys- 
tems, these unusual properties stem from the interplay 
between electronic, magnetic, lattice and orbital degrees 
of freedom. Understanding the microscopic process in- 
volved in these materials and the interaction of electrons 
with other entities is essential to understand the mecha- 
nism of high temperature superconductivity. 

Since the discovery of high-Tc superconductivity in 
cupratesflj, angle-resolved photoemission spectroscopy 
(ARPES) has provided key experimental insights in 
revealing the electronic structure of high temperature 
superconductors l)l3t ll4Hl5|) . These include, among oth- 
ers, the earliest identification of dispersion and a large 
Fermi surface lllq) . an anisotropic superconducting gap 
suggestive of a d— wave order parameter l|17|) . and an ob- 
servation of the pseudogap in underdoped samples ijlMl . 
In the mean time, this technique itself has experienced 
a dramatic improvement in its energy and momentum 
resolutions, leading to a series of new discoveries not 
thought possible only a decade ago. This revolution of 
the ARPES technique and its scientific impact result 
from dramatic advances in four essential components: 
instrumental resolution and efhciency, sample manipu- 
lation, high quality samples and well-matched scientific 
issues. 

The purpose of this treatise is to go through the 
prominent results obtained from ARPES on cuprate su- 
perconductors. Because there have been a number of 
recent reviews on the electronic structures of high-Tc 
materials ifla Il4t Hal , we will mainly present the latest 
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FIG. 1 Schematic of angle-resolved photoemission spec- 
troscopy. 



results not covered previously, with a special attention 
given on the electron-phonon interaction in cuprate su- 
perconductors. What has emerged is rich information 
about the anomalous electron-phonon interaction well 
beyond the traditional views of the subject. It exhibits 
strong doping, momentum and phonon symmetry de- 
pendence, and shows complex interplay with the strong 
electron-electron interaction in these materials. 



II. ANGLE-RESOLVED PHOTOEMISSION 
SPECTROSCOPY 



sample surface. Because of the symmetry breaking near 
the sample surface, the momentum perpendicular to the 
sample surface is not conserved. In order to obtain k^, 
one has to consider the inner potential which can be ob- 
tained in various ways (^19) . For strictly 2-dimensional 
materials or quasi-2-dimensional materials such as the 
cuprate superconductors discussed in this treatise, to the 
first approximation, one may treat k_L as a secondary 
effect. However, one should always be wary about the 
residual 3-dimensionality in these materials and its effect 
on photoemission data pOl) . 

The photoemission process can be understood intu- 
itively in terms of a "three step model" |21|): (i) Excita- 
tion of the electrons in the bulk by photons, (ii) Trans- 
port of the excited electrons to the surface, (iii) Emis- 
sion of the photoelectrons into vacuum. Under the "sud- 
den approximation" (described below), photoemission 
measures the single— particle spectral function A(k,tj), 
weighted by the matrix element M and Fermi function 
f(cj): I-A(k,w)|M|2f(w)|2l ^. The matrix element 
|Mp term indicates that, besides the energy and mo- 
mentum of the initial state and the final state, the mea- 
sured photoemission intensity is closely related to some 
experimental details, such as energy and polarization of 
incident light, measurement geometry and instrumental 
resolution. The inclusion of the Fermi function accounts 
for the fact that the direct photoemission measures only 
the occupied electronic states. 

The single-particle spectral function A(k,w) can be 
written in the following way using the Nambu-Gorkov 
formalism: 



A. Principle 

Angle-resolved photoemission spectroscopy is a pow- 
erful technique for studying the electronic structure of 
materials(Fig. [T|l 1)191) . The information of interest, i.e., 
the energy and momentum of electrons in the material, 
can be inferred from that of the photoemitted electrons. 
This conversion is made possible through two conserva- 
tion laws involved in the photoemission process: 
(1). Energy conservation: Es—hi'-Ekin-^', 
(2). Momentum conservation: K||=k||-|-G. 
where E^ represents the binding energy of electrons in 
the material; hi/ the photon energy of incident light; Fikin 
the kinetic energy of photemitted electrons; $ work func- 
tion; k|| momentum of electrons in the material paral- 
lel to sample surface; K|| projected component of mo- 
mentum of photoemitted electrons on the sample sur- 
face which can be calculated from the kinetic energy by 
hK\\=^/2rnEki^smd with h being Planck constant; G re- 
ciprocal lattice vector. Therefore, by measuring the in- 
tensity of the photoemitted electrons as a function of the 
kinetic energy at different emission angles, the electronic 
structure of the material under study, i.e., energy and 
momentum of electrons, can be probed direct ly lfl9l) . 

For 3-dimensional materials, the electronic structure 
also relies on kj^, the momentum perpendicular to the 



A{k,uj) = -{l/Tr)ImGii{k,uj) 



(2.1) 



G(fc,w) 



Z(fc, uj)ujTo + (e(fc) + x{k, uj))t2 + cj){k, uj)ti 



{Z(k, lj)ojY - is{k) + x(fc, u;)y - ^{k, toy 

(2.2) 
where Z, %, and (p represent a renormalization due to 
either electron-electron or electron-phonon interactions 
and e{k) is the bare-band energy, tq, ri, T2 are the matri- 
ces, and Gil represents the Pauli electronic charge den- 
sity channel measured in photoemission. In the weak 
coupling case, Z=l, x — 0, and (j) — A, the supercon- 
ducting gap. The same formalism can be extended to 
the normal state by setting cf> = 0. In the normal state, 
the spectral function can be written in a more compact 
wav ( , 22: 23) . in terms of the real and imaginary parts of 
the electron self energies ReE and ImE: 



A{k,uj) = - 



|ImI](fc,w)| 



tt(uj- e(k) - ReS(/c, uj))^ + (ImI](A:, w))2 

(2.3) 
where ReS describes the renormalization of the disper- 
sion and ImS describes the lifetime. 

In relating the photoemission process in terms of 
single particle spectral function A(k,a'), it is helpful to 



recognize some prominent assumptions involved: 

(1). The excited state of the sample (created by the 
ejection of the photo-electron) does not relax in the 
time it takes for the photo-electron to reach the detec- 
tor. This so-called "sudden-approximation" allows one to 
write the final state wave-function in a separable form, 
\]/^ = <i)^v[/^~i^ where $^ denotes the photoelectron and 

"^ f~^ denotes the final state of the material with N-1 
electrons. If the system is non-interacting, then the final 
state overlaps with a single eigenstate of the Hamiltonian 
describing the N-1 electrons, revealing the band structure 
of the single electron. In the interacting case, the final 
state can overlap with all possible eigenstates of the N-1 
system. 

(2) In the interacting case, A{k,uj) describes a "quasipar- 
ticle" picture in which the interactions of the electrons 
with lattice motions as well as other electrons can be 
treated as a perturbation to the bare band dispersion, 
e{k), in the form of a self energy, I](fc,u;). The validity 
of this picture as well as (1) rests on whether or not the 
spectra can be understood in terms of well-defined peaks 
representing poles in the spectral function. 
(3). The surface is treated no differently from the bulk 
in this A(A:,w). In reality surface states are expected 
and are observed and can lead to confusion in the data 
interpretation!!^. Surface termination also affects pho- 
toemission process |24). 

In addition to the matrix element M, there are other 
extrinsic effects which contribute to measured photoe- 
mission spectrum, e.g., the contribution from inelastic 
electron scattering. On the way to get out from inside 
the sample, the photoemitted electrons will experience 
scattering from other electrons, giving rise to a relatively 
smooth background in the photoemission spectrum. 



B. Technique 

As seen in Fig^ an ARPES system consists of a light 
source, chamber and sample manipulation and char- 
acterization systems, and an electron energy analyzer. 
Fig|2 is an example of a modern ARPES setup with the 
following primary components: 

(1). Light source: Possible light sources for angle- 
resolved photoemission are X-ray tubes, gas-discharge 
lamps, synchrotron radiation source and VUV lasers. 
Among them, the synchrotron radiation source is the 
most versatile in that it can provide photons with 
continuously tunable energy, fixed or variable photon 
polarization, high energy resolution and high photon 
fiux. The latest development of the VUV laser is 
significant as a result of its super-high energy resolution 
and super-high photon flux. In addition, the lower 
photon energy achievable by the VUV lasers makes the 
measured electronic structure more bulk-sensitive in 
certain materials (|25i|l . However, the strong final state ef- 




FIG. 2 A representative ARPES system on Beamline 10.0.1 
at the Advanced Light Source, Lawrence Berkeley National 
Lab. 



feet may limit its application to certain material systems. 

(2). Chambers and sample manipulation and char- 
acterization systems: In most of the photon energy 
range commonly used (20~100 eV), the escape depth 
of photoemitted j^^ctrons is on the order of 5~20 A, 
as seen in Fig|21:2G). This means that photoemission 
is a surface-sensitive technique. Therefore, obtaining 
and retaining a clean surface during measurement is 
essential to probe the intrinsic electronic properties of 
the sample. To achieve this, the ARPES measurement 
chamber has to be in ultra-high vacuum, typically better 
than 5x10^"'^^ Torr. A clean surface is usually obtained 
either by cleaving samples in situ in the chamber if the 
samples are cleavable or by sputtering and annealing 
process if the sample is hard to cleave. The quality of 
the surface can be characterized by low energy electron 
diffraction (LEED) or other techniques such as scanning 
tunneling microscopy (STM). The sample transfer 
system is responsible for quickly transferring samples 
from air to UHV chambers while not damaging the 
ultra-high vacuum. The manipulator is responsible for 
controlling the sample position and orientation, it also 
holds a cryostat that can change the sample temperature 
during the measurement. An advanced low temperature 
cryostat which can control the sample temperature 
precisely and has multiple degrees of translation and 
rotation freedoms is critical to an ARPES measuremnet. 

(3). Electron energy analyzer: An analyzer mea- 
sures the intensity of photoemitted electrons as a func- 
tion of their kinetic energy, i.e.. Energy Distribution 
Curve(EDC), at a given angle relative to the sample 
orientation. The dramatic improvement of the ARPES 



technique in the last decade is in large part due to the 
advent of modern electron energy analyzer, in particular, 
the Scienta series hemisphere analyzers. The enhance- 
ment of the performance lies in mainly three aspects: 
(i). Energy resolution improvement. 
The energy resolution of the electron energy analyzer 
improves steadily over time. The upgrade of the one- 
dimensional multichannel detection scheme of the VSW 
analyzer allows efficient measurement with ~20meV en- 
ergy resolution. Among others, it enabled the discovery 
of the d-wave superconducting gap structure l(l7l) . The 
first introduction of the Scienta 200 analyzer in the mid- 
dle 1990's dramatically improved the energy resolution 
to better than 5 meV. The latest Scienta R4000 analyzer 
has improved the energy resolution further to better than 
1 meV, as seen in Fig |^25') . 

We note that the total experimental energy resolution 
relies on both the analyzer resolution and the light source 
resolution. Sample temperature can also cause thermal 
broadening which is a limitation in some cases. The ne- 
cessity of multiple degrees of rotation controls as well as 
the exposure of the surface during an ARPES measure- 
ment often puts a lower limit on the sample tempera- 
ture. In addition, one should be aware of some intrinsic 
effects associated with the photoemission process, i.e., 
space charge effect and mirror charge effect (2r). When 
pulsed light is incident on a sample, the photoemitted 
electrons experience energy redistribution after escap- 
ing from the surface because of the Coulomb interaction 
between them (space charge effect) and between photo- 
emitted electrons and the distribution of mirror charges 
in the sample (mirror charge effect). These combined 
Coulomb interaction effects give rise to an energy shift 
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FIG. 3 Escape depth of photoemitted electron as a function 
of kinetic energy! 26) . For elements and inorganic compounds, 
the escape depth is found to follow the "universal curve" (red 
solid line). 
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FIG. 4 (a). Ultrahigh-resolution photoemission spectrum of 
an evaporated gold film measured using Scienta R4000 ana- 
lyzer at a temperature of 2.9 K (red circles), together with the 
Fermi-Dirac function at 2.9 K convolved by a Gaussian with 
full width at half maximum of 360 ^eV (a blue line). Total en- 
ergy resolution of 360 /ieV was confirmed from the very good 
match between the experimental and calculated spectral 25). 
The energy resolution from the VUV laser is estimated to be 
260 ^eV. (b). Angle mode testing image of Scienta R4000 
electron analyzer, the test was performed using "wire-and- 
slit" setup, with the angle interval between adjacent slits be- 
ing 1 degrees. In this particular angular mode, the analyzer 
collects emission angle within 30 degrees simultaneously. 



and a broadening whose magnitude depends on the pho- 
ton energy, photon flux, beam spot size, emission angles 
and etc. For a typical third-generation synchrotron light 
source, the energy shift and broadening can be on the 
order of 10 meV (FigEJllj). This value is comparable to 
many fundamental physical parameters actively studied 
by photoemission spectroscopy and should be taken seri- 
ously in interpreting photoemission data and in designing 
next generation experiments. 

(ii). Momentum resolution; 
The introduction of the angular mode operation in the 
new Scienta analyzers has also greatly improved the an- 
gular resolution, from a previous ~2 degrees to 0.1~0.3 
degree. This improvement of the momentum resolu- 
tion allows one to observe detailed structures in the 
band structure and Fermi surface, as well as subtle 
but important many-body effects. As an example, re- 
cent identification of two Fermi surface sheets (so-called 
"bilayer splitting" )in BiaSraCaCuaOg (Bi2212) (Fig© 
is largely due to such an improvement of momentum 
resolution((2a |2a l^fl) , combined with the advancement 
of theoretical calculations l|24J) . 

(iii). Two-dimensional multiple angle detection; 
Traditionally, the electron energy analyzer collects one 
photoemission spectrum, i.e., energy distribution curve 
(EDC), at one measurement for each emission angle. 
Modern electron energy analyzers collect multiple angles 
simultaneously. As shown in Fig^, the latest Scienta 
R4000 analyzer can collect photoemitted electrons in the 
angle range of 30 degrees simultaneously. Therefore, at 
one measurement, the raw data thus obtained, shown 
in Fig[7^, is a 2-dimensional image of the photoelectron 
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FIG. 5 Space charge and mirror charge effect in 
photoemission i27r) . Fermi edge broadening (solid square) and 
the Fermi edge shift (open circle) as a function of sample cur- 
rent. The beam spot size is ~0.43mmx0.30 mm. The inset 
shows the measured overall Fermi edge width as a function of 
the sample current, which includes all contributions including 
the beamline, the analyzer and the temperature broadening. 
The net broadening resulting from pulsed photons is obtained 
by deconvolution of the measured data, taking the width at 
low photon flux as from all the other contributions. 



tribution curves (EDCs), which represent photoelectron 
intensity as a function of energy for a given momentum. 
The 2D image comprising the raw data is then equivalent 
to a number of EDCs at different momenta (FigjTls). The 
peak position at different momenta will give the energy- 
momentum dispersion relation determining the real part 
of electron self-energy ReE. The EDC linewidth deter- 
mines the quasiparticle lifetime, or the imaginary part of 
electron self-energy ImE. However, the EDC lineshape is 
usually complicated by a background at higher binding 
energy, the Fermi function cutoff near the Fermi level, 
and an undetermined bare band energy which make it 
difhcult to extract the electron self-energy precisely. 



An alternative way to visualize the 2D data is to ana- 
lyze photoelectron intensity as a function of momentum 
for a given electron kinetic energy (|33)by means of mo- 
mentum distribution curves (MDCs) (J33 : yj). This ap- 
proach provides a different way of extracting the elec- 
tron self-energy. As shown in Fig[7|:, the MDCs exhibit 
well-defined peaks with flat backgrounds; moreover, they 
can be fitted by a Lorentzian lineshape. When the band- 
width is large, the band dispersion ek can be approx- 
imated as ek = vok in the vicinity of the Fermi level. 
Under the condition that the electron self-energy shows 
weak momentum dependence, A(k,Ci;) indeed exhibits a 
Lorentzian lineshape as a function of k for a given binding 
energy . By fitting a series of MDCs at different binding 
energies to obtain the MDC position k and width F (full- 
width at half maximum, FWHM) (Fig[7|i)@): one can 
extract the electron self-energy directly as: HeTi—huj-kYo 
and ImI]=Fvo/2. 
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FIG. D (a.). Jixper 
Pb-doped Bi2212J33) 
612212(12^ 



intensity (represented by false color) as a function elec- 
tron kinetic energy and emission angle (and hence mo- 
mentum). This 2-dimensionality greatly enhances data 
collection efficiency and provides a convenient way of an- 
alyzing the photoemission data. 

As shown in FigEI , the traditional way to visualize the 
photoemission data is by means of so-called energy dis- 



It is worthwhile to point out the latest effort in at- 
tempting to overcome the surface sensitivity issue related 
with photoemission. As seen from Fig|31 in the usual 
photon energy range used for valence band photoemis- 
sion, the photoemitted electron escape depth is on the 
order on 5^10 A. Therefore, it is always an issue whether 
the photoemission results obtained in this energy range 
represents the bulk properties. To overcome such a prob- 
lem, there have been two approaches by employing either 
high photon energy or lower photon energy. As seen from 
Figl^l when the photon energy is on the order of 1 KeV, 
the electron escape depth can be increased to ^^20 A(36). 
However, this modest enhancement of the bulk sensitiv- 
ity comes at a price of sacrificing both the energy res- 
olution and momentum resolution. On the other hand, 
when the photon energy is low, one can see that the elec- 
tron escape depth increases dramatically. Note that this 
"universal" curve is obtained from metals, whether the 
same curve can be applied to oxide materials remains un- 
clear yet. In addition to the potential engancement of the 
bulk sensitivity, one may further improve the energy and 
momentum resolution by going to lower photon energy.. 
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FIG. 7 Illustration of the MDC method for extracting the 
electron self-energy, (a) Raw photoemission data for LSCO 
with x=0.063 (Tc--12K) along the (0,0)-(7r,7r) nodal direction 
at 20 Kya). The two-dimensional data represent the photo- 
electron intensity (denoted by false color) as a function of en- 
ergy and momentum, (b) Energy distribution curves (EDCs) 
at different momenta. The EDC colored red corresponds to 
the Fermi momentum kp- (c) Momentum distribution curves 
(MDCs) at different binding energies. The MDC colored red 
corresponds to the Fermi level, (d) Energy-momentum dis- 
persion relation extracted by the MDC method. The inset 
shows the MDC width as a function of energy. 



III. ELECTRONIC STRUCTURES OF HIGH 
TEMPERATURE SUPERCONDUCTORS 





FIG. 8 (a) Schematic of the real-space Cu02 plane. The 
Cu02 plane consists of copper (pink solid circles) and oxygen 
(black open circles), (b) The corresponding Brillouin zone 
in a reciprocal space. In the first Brillouin zone, the area 
near {n/2, 7r/2) (denoted as red circle) is referred to as nodal 
region, and the (0,0)-(7r,7r) direction is the nodal direction (red 
arrow). The area near (■7r,0) and (0,7r) is referred to as the 
antinodal region (shaded circles) . The blue solid line shows a 
schematic Fermi surface. 




FIG. 9 Bonding in Cu02 plane(|4Q). The atomic Cu 3d level 
is split due to the cubic crystal field into e^ and t25 states. 
There is a further splitting due to an octahedral crystal field 
into x^ - y^, 3z^ - r^, xy, and xz, yz states. For divalent Cu 
which has nine 3d electrons, the uppermost x^ - y^ level is 
half filled, while all other levels are completely filled. There is 
a strong hybridization of the Cu states, particularly the x^ - 



A. Basic Crystal Structure and Electronic Structure 

A common structural feature of all cuprate supercon- 
ductors is the Cu02 plane (FiglSJi) which is responsible 
for the low lying electronic structure; the Cu02 planes 
are sandwiched between various block layers which serve 
as charge reservoirs to dope Cu02 planes l)37t l38D . For 
the undoped parent compound, such as La2Cu04, the 
valence of Cu is 2+, corresponding to 3d^ electronic con- 
figuration. Since the Cu^^ is surrounded by four oxygens 
in the Cu02 plane and apical oxygen(s)or halogen(s) per- 
pendicular to the plane, the crystal field splits the oth- 
erwise degenerate five rf-orbitals, as schematically shown 
in Fig|ni39). The four lower energy orbitals, including 
xy, xz, yz and 3z^ - r^, are fully occupied, while the or- 
bital with highest energy, x^-y^, is half-filled, since the 
energies of the Cu d-orbitals and O 2p-orbitals are close, 
there is a strong hybridization between them. As a re- 
sult, the topmost energy level has both Cu (1^2 _y2 and O 



dimensional Cu 3d^2_y2-0 2pi,y antibonding dpo" band. The 
hybridization of the other 3d levels is smaller and is indicated 
in Figure only by a broadening. 



2px,y character. 

The same conclusion is also drawn from band struc- 
ture calculations (Figll0b)l39l). According to both 
simple valence counting (Fig|3 and band structure 
calculation (Fig llUb ). the undoped parent compound 
is supposed to be a metal. However, strong Coulomb 
interactions between electrons on the same Cu site 
makes it an antifcrromagnctic insulator with an energy 
gap of 2 eV(42; 43). The basic theoretical model for the 
electronic structure most relevant to our discussion is 
the multi-band Hubbard Hamiltonian l|44i |45|) containing 
d states on Cu sites, p states on O sites, hybridiza- 
tion between Cu-0 states, hybridization between 0-0 
states, and Coulomb repulsion terms. In terms of hole 



notation, i.e., starting from the filled-shell configuration 
(3d^*',2p^) corresponding to a formal valence of Cu^+ and 
O^^, the general form of the model can be written as(J4y): 



H = ^ Sdd-l^dia- + ^ SpviaPia + ^ tpdPtjdi^ + h. 



la 



<li>(j 



^ UdUi^nii + ^ to-optaPi'<y + h-^- 

i <ll'>a 



y^ UpdnicrUicy' + ^ UpTii^nii 

<il>(T(T' I 



(3.1) 



where the operator d^ creates Cu (3d2.2_j,2)holcs at site 
i, and p^ creates 0(2p) holes at the site I. Ud is the on- 
site Coulomb repulsion between two holes on a Cu site. 
The third term accounts for the direct overlap between 
Cu-0 orbitals. The fifth terms describes direct hopping 
between nearest-neighbor oxygens, and \Jpd in the sixth 
term is the nearest neighbor Coulomb repulsion between 
holes on Cu and O atoms. Qualitatively, this model gives 
the energy diagram in FigllOb. 

Simplified versions of model Hamiltonians have also 
been proposed. Notably among them are the single-band 
Hubbard model(47) and t-J model(^). The t-J Hamil- 
tonian can be written in the following forml)4a: l50|) : 



Htj^ 



/ , y'^ia'^jo 



-t > (5„ 

<ij>:0- 



Cj(j -\- 



iJ.c.) + J ^(Si.Sj-n,Tn,i/4) 



<i]> 



;^t 



.t 



(3.2) 
where the operator cl^=cl^{l — fii-a) excludes double 
occupancy, J = 4t^/[/ is the antiferromagnetic exchange 
coupling constant, and S^ is the spin operator. Since 
the hopping process may also involve the second (t ) and 
third (t ) nearest neighbor, an extended t-J model, the 
t— t — t"— J model, has also been proposed ((sj). 



B. Brief Summary of Some Latest ARPES Results 

ARPES has provided key information on the electronic 
structure of high temperature superconductors, including 
the band structure, Fermi surface, superconducting gap, 
and pseudogap. These topics are well covered in recent 
reviewsfl4"; 15) that we will not repeat here. Instead, we 
briefly summarize some of the latest developments not 
included before. 

Band structure and Fermi surface: The bi-layer split- 
ting of the Fermi surface is well established in the over- 
doped Bi2212(J2i m 113), as shown in FigEl and also 
suggested to exist in underdoped and optimally doped 
Bi2212(|53 ISa IM: ISq) . Recent measurements also show 
that there is a slight splitting along the (0,0)-(7r,7r) 
nodal direction l|56|) . The measurement on four-layered 
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FIG. 10 (a). LDA calculated band structure of La2CuO4(0). 



The band labeled B is bonding band between Cu 3d^ 



and 



O 2p states while the band labeled A is the corresponding 
antibonding band that is half-filled; (b). Schematic of Zhang- 
Rice singlet statedia: 143) . (c). Schematic energy diagrams 
for undoped and doped Cu02 planes (1421 . (cl). Band picture 
for a half-filled (undoped) Cu02 plane (Fermi liquid); (c2). 
Charge-transfer insulating state of the Cu02 plane with split 
Cu 3d bands due to on-site Coulomb repulsive interaction U. 
The O 2p band is separated by a charge transfer energy A 
from the upper Cu 3d band; (c3) and (c4) show rigid charge 
transfer energy bands doped with holes and electrons, respec- 
tively; (c4). Formation of mid-gap states inside the charge 
transfer gap. 



Ba2Ca3Cu408F2 has identified at least two clear Fermi 
surface sheets ((5 7^). 

Superconducting gap and pseudogap: Since the first 
identification of an anisotropic superconducting gap in 
Bi2212p7f). subsequent measurements on the supercon- 
ductors such as Bi2212||5i|52|63;j6D), Bi2201(62; 63), 
Bi2223(|6l m M), YBaaCuaOT-alEl, LSCO(68) have 
established a universal behavior of the anisotropic su- 
perconducting gap in these hole-doped superconductors 
which is consistent with d-wave pairing symmetry (al- 
though it is still an open question whether the gap form 
is a simple d-wave-like A(/c)=Ao[cos(k2;a)-cos(kj^a)] or 
higher harmonics of the expansion should be included) . 
The measurements on electron-doped superconductors 
also reveal an anisotropic superconducting gaD l)69H70|) . 

One interesting issue is, if a material has multiple 
Fermi surface sheets, whether the superconducting gap 
on different Fermi surface sheets is the same. This is- 
sue traces back to superconducting SrTiOa where it was 
shown from tunneling measurements that different Fermi 
surface sheets may show different Fermi surface gapsfTr). 
With the dramatic advancement of the ARPES tech- 
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nique, different superconducting gaps on different Fermi 
surface sheets have been observed in 2H-NbSeT 1172) and 
MgB2(,73). For high-T^ materials, Bi2212 shows two 
clear FS sheets, but no obvious difference of the supercon- 
ducting gas has been resolved |6l|) . In Ba2Ca3Cu408F2, 
it has been clearly observed that the two Fermi surface 
sheets have different superconducting gapsfS^)- 

Time reversal symmetry breaking: It has been pro- 
posed theoretically that, by utilizing circularly polarized 
light for ARPES, it is possible to probe time-reversal 
symmetry breaking that may be associated with the 
pseudogap state in the underdoped samples IJT^ iTSf) . 
Kaminski et al. first reported the observation of such 
an effect l(76f l. However, this observation is not repro- 
duced by another group l)77j) and the subject remain 
controversial ifrSjl . 



IV. ELECTRON-PHONON COUPLING IN HIGH 
TEMPERATURE SUPERCONDUCTORS 

The many-body effect refers to interactions of electrons 
with other entities, such as other electrons, or collective 
excitations like phonons, magnons, and so on. It has 
been recognized from the very beginning that many-body 
effects are key to understanding cuprate physics. Due 
to its proximity to the antiferromagnetic Mott insulat- 
ing state, electron-electron interactions are extensively 
discussed in the literature l(l4L llSJ) . In this treatise, we 
will mostly review the recent progress in our understand- 
ing of electrons interacting with bosonic modes, such 
as phonons. This progress stems from improved sam- 
ple quality, instrumental resolution, as well as theoreti- 
cal development. In a complex system like the cuprates, 
it is not possible to isolate various degrees of freedom 
as the interactions mix them together. We will dis- 
cuss the electron-boson interactions in this spirit, and 
will comment on the interplay between electron-phonon 
and electron-electron interactions whenever appropriate. 
Here by bosonic modes, we are referring to collective 
modes with sharp collective energy scale such as the op- 
tical phonons and the famous magnetic resonance mode 
seen in some cupratesJjO; 80; 81), but not the broad exci- 
tation spectra such as those from the broad electron/spin 
excitations as these issues have been discussed in previous 
reviews. Furthermore, we believe the effects due to sharp 
mode coupling seen in cuprates are caused by phonons 
rather than the magnetic resonance. Our reason for not 
attributing the observed effect to magnetic resonance will 
become apparent from the rest of the manuscript. With 
more limited data, other groups have taken the view that 
the magnetic resonance is the origin of the boson cou- 
pling effect. For this reason, we will focus more on our 
own results in reviewing the issues of electron-phonon 
interaction in cuprates. 

The electron-phonon interactions can be characterized 
into two categories: (i). Weak coupling where one can 
still use the perturbative self-energy approach to describe 



the quasiparticle and its lifetime and mass; (ii). Strong 
coupling and polaron regime where this picture breaks 
down. 



A. Brief Survey of Electron-Phonon Coupling in 
High-Temperature Superconductors 

It is well-known that, in conventional superconductors, 
electron-phonon (cl-ph) coupling is responsible for the 
formation of Cooper pairs(4). The discovery of high 
temperature superconductivity in cuprates was actually 
inspired by possible strong electron-phonon interaction 
in oxides owing to polaron formation or in mixed- valence 
systems (|l|). However, shortly after the discovery, a 
number of experiments lead some people to believe that 
electron-phonon coupling may not be relevant to high 
temperature superconductivity. Among them are l)82() : 
(1). High critical transition temperature T^ 

So far, the highest Tc achieved is 135 K in 
HgBa9Ca9Cu^O« l|83|) at ambient pressure and ^160 
K under high pressure |84l) . Such a high Tc was not 
expected in simple materials using the strongly coupled 
version of BCS theory, or the McMillan equations. 
(2). Small isotope effect on Tc 

It was found that the isotope effect in optimally- 
doped samples is rather small, much less than 
that expected for strongly-coupled phonon-mediated 
superconductivity (|85|) . 
(3). Transport measurement 

The linear resistivity-temperature dependence in 
optimally doped samples and the lack of a saturation 
in resistivity over a wide temperature range have been 
taken as an evidence of weak electron-phonon coupling 
in the cuprate superconductors ijSql . 
(4) . d-wave symmetry of the superconducting gap 

It is generally believed that electron-phonon coupling 
is favorable to s-wave coupling. 
(5). Structural instability. 

It is generally believed that sufficiently strong 
electron-phonon coupling to yield high Tc will result in 
structural instabilityljsj) . 

Although none of these observations can decisively rule 
out the electron-phonon coupling mechanism in high- 
Tc superconductors, overall they suggest looking else- 
where. Instead, strong electron-electron correlation has 
been proposed to be the mechanism of high-Tc supercon- 
ductivity 1|88|) . This approach is attractive since d— wave 
pairing is a natural consequence. Furthermore, the high 
temperature superconductors evolve from antiferromag- 
netic insulating compounds where the electron-electron 
interactions are strong (Q; |3) 

However, there is a large body of experimental evi- 
dence also showing strong electron-phonon coupling in 
high-temperature superconductors l|89t lOOt Wof . Among 
them are: 
(1). Isotope effect; 
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FIG. 11 Doping dependence of the oxygen isotope effect ao 
on Tc in several classes of cuprates i92t 193 : [941) . The "1/8 
anomaly" data found in LSCO system is highlighted in the 
shaded region. 
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FIG. 12 Schematic of Big mode (a), Aig mode (b), half- 
breathing mode (c), full-breathing mode (d) and apical oxy- 
gen mode (e). 



As seen in Fig llll although at the optimal doping, 
the oxygen isotope effect on Tc is indeed small, it get s 
larger and becomes significant with reduced doping((93). 
In particular, near the "1/8" doping level, the isotope 
effect in (La2-2:Sra;)Cu04 and (La2-xBa2;)Cu04 is 
anomalously strong, which is related to the structural 
inst ability l|94|) . Furthermore, the measurement of an 
oxygen isotope effect on the in-plane penetration depth 
also suggests the importance of lattice vibration for 
high-Tc superconductivtvl|95r). 

(2). Optical spectroscopy and Raman scattering; 

Raman scattering 119^ and infrared spectroscopv l|97J) 
reveal strong electron-phonon interaction for certain 
phonon modes. Some typical vibrations related to the 
in-plane and apical oxygens are depicted in Fig ll2l In 
YBa2Cu307_5, it has been found that, the Big phonon, 



FIG. 13 Anomalou s soft ening of the Big phonon when YBCO 
is cooled below T^ ilOOD . The inset shows the fit of a Fano 
function to the phonon peak at T=72K l98r) 



which is related to the out-of-plane, out-of-phase, in- 
plane oxygen vibrations (see Fig ll2|l . exhibits a Fano-like 
lineshape (Fig^J and shows an abrupt sof tenin g upon 
entering the superconducting state l|98i lOOl : llOOf) . The 
Aig modes, as found in HgBa2Ca.ciC u40in (Hgl234)(|l0l|) 
and in HgBa2Ca2Cu308 ('Hgl223) |l0^ . exhibit espe- 
cially strong superconductivity-induced phonon soften- 
ingfFig [n|) . Infrared reflectance measurements on var- 
ious cuprates found that the frequency of the Cu-0 
stretching mode in the Cu02 plane is v ery sensitive to 
the distance between copper and oxvgen |97|) . 

Fig. El shows Raman data as a function of doping 
in LSCO 1103). The sharp structures at high frequency 
are signals from multiphonon processes, which can only 
occur if the electron-phonon interaction is very strong. 
One can see that this effect is very strong in undoped and 
deeply underdoped regime, and gets weaker with doping 



increase. 



(3). Neutron scattering 

Neutron scattering measurements have provided rich 
information about electron-phonon coupling in high tem- 
perature superconductors! 104; 105; 106). As seen from 
FigllBb. the in-plane "half-breathing" mode exhibits 
strong frequency renormalizations upon doping along 
(001) direction! 104; 107). In (Lai.85Sro.i5)Cu04, it is re- 
ported that, at low temperature, the half-breathing mode 
shows a discontinuity in dispersion fFigllGbldOSl. In 
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YBCO, neutron scattering indicates that the softening of 
the Big mode upon entering the superconducting state 
is not just restricted near q=0, as indicated by Raman 
scattering fFig ll3|l . but can be ob served in a large part 
of the BriUouin zone fFig [r7Jl|T06|) . 

(4). Material and structural dependence; 

There is a strong material and structural depen- 
dence to t he h i gh-T r superconductivity, as exemplified 
in Figtra ifiol lUOl) . Empirically it is found that, for 
a given homologous series of materials, the optimal Tc 
varies with the number of adjacent Cu02 planes, n, in a 
unit cell: Tc goes up first with n, reaching a maximum 
at n=3, and goes down as n further increases. For the 
cuprates with the same number of Cu02 layers, Tc also 
varies significantly among different classes. For exam- 
ple, the optimal Tc for one-layered (La2-a;Sra;)Cu04 is 
40K while it is 95K for one-layered HgBa2Cu04. These 
behaviors are clearly beyond simplified models that con- 
sider Cu02 planes only, such as the t— J model. In fact, 
such effects were taken as evidence against theoretical 
models based on such sim ple m odels and in favor of the 
interlayer tunneling model(lil2). Although the interlayer 
tunneling model has inconsistencies with some experi- 
ments, the issue that the material dependence cannot be 
explained by single band Hubbard and t-J model remains 
to be true. 

The above results suggest that the lattice degree of 
freedom plays an essential role. However, the role of 
phonons has not been scrutinized as much, in particu- 
lar in regard to the intriguing question of whether high- 
Tc superconductivity involves a special type of clcctron- 
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FIG. 14 Raman spectra of Hgl234 showing a giant 
superconductivity-induced mode softening across Tc=123 
K ( 101 ). The modes at 240 cm~^ and 390 cm~^ correspond to 
Aig out-of-plane, in-phase vibration of oxygens in the Cu02 
planes. Upon cooling from room temperature to 4.5 K, the 
240 cm~^ Aig mode shows a abrupt drop in frequency at Tc 
from 253 to 237 cm~ ^ and the 390 cm"'^ mode drops from 
395 to317cm-^((Tol. 
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FIG. 15 Aig two-phonon Raman spectra in LSCO at different 
dopings. The dark gray area indicates that the two-phonon 
peak of the (7r,7r) LO mode is strong and the light gray area 
indica tes th at the two-phonon peak of the ( 7r,0) LO mode is 
strongJlQ3)- 
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FIG. 16 (a). Dispersion of the Cu-O bond-stret ching vibra- 
tions in the (lOO)-direction in (La2-a;Sra;)Cu04(|io3)- (b). 
Anomalous dispersion of LO phonons in Lai.85Sro.i5Cu04. 
lOK data are filled circles and room temperature data are 
empty squares. Grey shaded circles indicate the frequency of 
the weak extra branch seen at 10K(jlO^. 
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FIG. 17 q dependence of Big mode peak position at differ- 
ent temperatures in YBCO. Dashed lines are guides to the 
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phonon coupling. In other words, the complexity of 
electron-phonon interaction has not been as carefully ex- 
amined as some of the electronic models. As a result, 
many naive arguments are used to argue against electron- 
phonon coupling as if the conclusions based on simple 
metals are applicable here. Recently, a large body of 
experimental results from angle-resolved photoemission, 
as we review below, suggest that electron-phonon cou- 
pling in cuprates is not only strong but shows behaviors 
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distinct from conventional electron-phonon coupling. In 
particular, the momemtum dependence and the interac- 
tion between electron-phonon interaction and electron- 
electron interaction are very important. 



B. Electron-Phonon Coupling: Theory 

1. General 

Theory of electron-phonon interaction in the presence 
of strong electron correlation has not been developed. 
Given both interactions are important in cuprates, it is 
difficult a priori to have a good way to address these 
issues. In fact, we believe that an important outcome of 
our research is the stimulus to develop such a theory. In 
the mean time, our strategy to separate the problem in 
different regimes and see to what extent we can develop 
a heuristic understanding of the experimental data. Such 
empirical findings can serve as a guide for comprehensive 
theory. We now start our discussion with an overview of 
existing theories of electron-phonon physics. 

The theories of electron-phonon coupling in condensed 
matter have been developed rather separately for met- 
als and insulators. In the former case, the dominant 
energy scale is the kinetic energy or the Fermi energy 
ep on order of 1 — lOeV, and the phonon frequency 
57 ^ 1 — lOOmeV is much smaller. The Fermi degeneracy 
protects the many-body fermion system from perturba- 
tions and only the small energy window near the Fermi 
surface responds. Therefore even if the lattice relaxation 
energy Elpi = g^ /lo for the localized electron is compa- 
rable to the kinetic energy ep the el-ph coupling is essen- 
tially weak and the perturbative treatment is justified. 
The dimensionless coupling constant A is basically the 
ratio of Elji/ef, which ranges A = 0.1 — 2 in the usual 
metals. In the diagrammatic language, the physics de- 
scribed above is formulated within the framework of the 
Fermi liquid theory (112 ). The el-el interaction is taken 
care of by the formation of the quasi-particle, which is 
well-defined near the Fermi surface, and the el-ph vertex 
correction is shown to be smaller by the factor of fl/ep 
and can be neglected. Therefore the multi-phonon exci- 
tations are reduced and the single-loop approximation or 
at most the self-consistent Born approximation is enough 
to capture the physics well, i.e., Migdal-Eliashberg for- 
malism. 

When a carrier is put into an insulator, on the other 
hand, it stays near the bottom of the quadratic disper- 
sion and its velocity is very small. The kinetic energy 
is much smaller than the phonon energy, and the carrier 
can be dressed by a thick phonon cloud and its effective 
mass can be very large. This is called the phonon po- 
laron. Historically the single carrier problem coupled to 
the optical phonon through the long range Coulomb in- 
teraction, i.e., Frohlich polaron, is the first studied model, 
which is defined in the continuum. When one considers 
the tight- binding models, which is more relevant to the 
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Bloch electron, the bandwidth W plays the role of ep 
in the above metallic case. Then again we have three 
energy scales, W, Erl, and 51. Compared with the 
metallic case, the dominance of the kinetic energy is not 
trivial, and the competition between the itinerancy and 
the localization is the key issue in the polaron problem, 
which is controlled by the diniensionless coupling con- 
stant A = Ejil/W . Another dimensionless coupling con- 
stant is 5 = Eul/H, which counts the number of phonon 
quanta in the phonon cloud around the localized electron. 
This appears in the overlap integral of the two phonon 
wavefunctions with and without the phonon cloud as: 



< phonon vacuum | phonon cloud >cx e 



(4.1) 



This factor appears in the weight of the zero-phonon line 
of the spectral function of the localized electron, and S 
can be regarded as the maximum value for the number 
of phonons A^ph near the electron. In a generic situa- 
tion, A'ph is controlled by A, and there are cases where 
iVph shows an (almost) discontinuous change from the 
itinerant undressed large polaron to the heavily dressed 
small polaron as A increases. This is called the self- 
trapping transition. Here a remark on the terminology 
"self-trapping" is in order. Even for the heavy mass po- 
laron, the ground state is the extended Bloch state over 
the whole sample and there is no localization. However 
a small amount of disorder can cause the localization. 
Therefore in the usual situation, the formation of the 
small polaron implies the self-trapping, and we use this 
language to represent the formation of the thick phonon 
clouds and huge mass enhancement. In cuprates, it is 
still a mystery why the transport properties of the heav- 
ily underdoped samples do not show the strong localiza- 
tion behavior even though the ARPES shows the small 
polaron formation as will be discussed in D.l. 

Now the most serious question is what is the picture for 
the el-ph coupling in cuprates ? The answer seems not 
so simple, and depends both on the hole doping concen- 
tration, momentum and energy. The half-filled undoped 
cuprate is a Mott insulator with antiferromagnetic order- 
ing, and a single hole doped into it can be regarded as the 
polaron subjected to the hole-magnon and hole-phonon 
interactions. At finite doping, but still in the antiferro- 
magnetic (AF) order, the small hole pockets are formed 
and the hole kinetic energy can be still smaller than the 
phonon energy. In this case the polaron picture still per- 
sists. The main issue is to what range this continues. 
One scenario is that once the antiferromagnetic order 
disappears the metallic Fermi surface is formed and the 
system enters the Migdal-Eliashberg regime. However, 
there are several physical quantities such as the resistiv- 
ity. Hall constant, optical conductivity, which strongly 
suggest that the physics still bears a strong characteris- 
tics of doped holes in an insulator rather than a simple 
metal with large Fermi surface. Therefore the crossover 
hole concentration x^ between the polaron picture and 
the Migdal-Eliashberg picture remains an open issue. 
Probably, it depends on the momentum/energy of the 



spectrum. For example, the electrons have smaller veloc- 
ity and are more strongly coupled to the phonons in the 
anti-nodal region near (±7r,0), (0,±7r), remaining pola- 
ronic up to higher doping, while in the nodal region, the 
electrons behave more like the conventional metallic ones 
since the velocity is large along this direction. Further- 
more, the low energy states near the Fermi energy are 
well described by Landau's quasi-particle and Migdal- 
Eliashberg theory, while the higher energy states do not 
change much with doping even at a; = 0.1 (|ll3) sugges- 
tive of polaronic behavior. In any event, the dichotomy 
between the hole doping picture and the metallic (large) 
Fermi surface picture is the key issue in the research of 
high Tc superconductors. 



2. Weak Coupling - Perturbative and Self-Energy Description 

We review first the Migdal-Eliashberg regime, in which 
the electron-phonon interaction results in single-phonon 
excitations and can be considered as a perturbation to 
the bare band dispersion. In this case, dominant features 
of the mode coupling behavior can be captured using the 
following form for the self-energy: 



Yj{k,U!) = T/N y^ g^{k, q)D{q, iiy)T3G{k - q.iui - iv)T-i 

q,u 

(4.2) 
-i?yi is the phonon propagator, 51^ is 

the phonon energy, T is temperature, N is the number of 
particles and t^ is the Pauli matrix. 

In this form of the self-energy, corrections to the 
electr on-ph onon vertex, g, are neglected as mentioned 
above ljllSJ) . Furthermore, we assume only one-iteration 
of the coupled self-energy and Green's function equa- 
tions. In other words, in the equation for the self- energy, 
E, we assume bare electron and phonon propagators. Go 
and Dq. With these assumptions, the imaginary parts of 
the functions Z, x, and 0, denoted as Z2, X2, and 02, 
are: 



where D{q,uj) 



2n„ 



Z2{k,u;)Lu = ^52(fc,g)(V2){[<5(c^ - n, - E^-g) 



+6{uj -ng + Ek-g)] [f(ng -Lo) + n{ng)] 



+[-s{Lj+ng-Ek-g)-s{Lo+ng+Ek-g)][fing+Lj)+ning)]} 

(4.3) 



X2ik,uj) = ^g^ik,q){Trek-q/2Ek-g){[-Siuj-ng-Ek-g) 



+ Siu ~ng+ Ek-gWiftg -L0)+ n{ng)] 
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FIG. 19 Self-energy for electrons coupled to an Einstein mode 
wit h O = 35 meV and electron-phonon vertex g = 0.15 
eVini). (al), (bl), and (cl) plots ImE= -Z2L0 + X2 for a 
normal state electron at 10 K, for a normal state electron 
at 100 K, and for an electron in an s-wave superconducting 
state at 10 K, respectively. (a2), (b2), and (c2) plots the 
corresponding real parts, ReE , obtained using the Kramers- 
Kronig relation. 
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hn{nq]} 
(4.5) 

where f{x), n{x), are the Fermi, Bosc distribution func- 
tions and Ek is the superconducting state dispersion. 
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The above equations are essentially those of Eliashberg 
theory for strongly-coupled superconductors. Although 
A can be large (> 1), i.e., "strongly-coupled", the vertex 
corrections and niulti-phonon processes are still negligi- 
ble due to the Fermi degeneracy and small 51/i?i?(116). 
To illustrate the essential features of mode coupling, 
we consider an Einstein phonon coupled isotropically 
to a parabolic band. We present this calculation in the 
spirit of Engelsberg and Schrieffer, who first calculated 
the spectral function for an electron-phonon coupled 



svstem l)ll7j) and which provided the foundation fo r the 
later work by Scalapino, Schrieffer, and Wilkins l)ll8|) 
in the superconducting state. Fig^l plots —Ziu: -I- X2, 
the imaginary part of the phonon self-energy, ImS, that 
represents the renormalization to the diagonal channel 
of the electron propagator, or the one in which the 
charge number density is subjected to electron-phonon 
interactions. This part of the self-energy gives a finite 
lifetime to the electron, and consequently broadens the 
peak in the spectra (ImE in Aik^uS) (Ea. l2.3|l is the 
half-width-at-half-maximum, HWHM of the peak). In 
the normal state, — Z2CJ + X2 takes the familiar form: 



Imi;(fc,w) = I]5-7rg^(fc,g)[2n(OJ-h 

/(r!,+cc;) + /(f7g-cc;)]<5(c^-£;fe_,) (4.6) 

which when integrated over q becomes: 

ImE(fc,tj) = / drJafcF(f7)[2n(r2) + /(f7 + w) + /(r2-w)] 

(4.7) 
where a^^,F(^)^ the Eliashberg function, represents the 
coupling of the electron with Fermi surface momentum 
fc, to all rj phonons connecting that electron to other 
points on the Fermi surface. 

For the normal state electron at lOK (Figll9bl). there 
is a sharp onset of the self-energy that broadens the spec- 
tra beyond the mode energy; for the normal state electron 
at lOOK (Fig. I19bl). the onset of the self-energy is much 
smoother and occurs over ~ 50meV; for the supercon- 
ducting state electron (Fig. I19bl). there is a singularity 
that causes a much more abrupt broadening of the spec- 
tra at the energy fi -I- A. The superconducting state 
singularity is due to the density of states pile-up at the 
gap energy; the energy at which the decay onsets shift 
by A, since below the gap energy there are no states to 
which a hole created by photoemission can decay. For 
each of these imaginary parts of the self-energy, one can 
use the Kramers-Kronig transform to obtain the real part 
of the self energy, which renormalizes the peak position 
(ReE in A(k, uS) (Ea. l2.3|l changes the position of the peak 
in the spectral function). The real self energies thereby 
obtained are also plotted in Fig. I19b2. Fig. I19b2. and 
Fig. [TOb2. In the superconducting state, again there is 
a singularity that causes a more abrupt break from the 
bare-band dispersion at the energy SI -I- A. 

For most metals, where the electrons are weakly inter- 
acting, and therefore the poles of the spectral function 
are well-defined, one would expect such a treatment to 
hold and indeed it does , as e videnced by several case s in- 
cluding Beryllium ini IT20|) and Molvbdenum lfm)) . A 
priori, one might not expect the same to hold in ce- 
ramic materials such as the copper-oxides, where the 
copper d-wave electrons are localized and subject to 
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FIG. 20 Simulated electron-phonon coupling using Einstein model. Spectral function (al, bl), EDCs (a2, b2), MDC-derived 
dispersion (a3, b3), and the MDC-derived width (a4, b4) (imaginary part of self-energy) for two different couplings (a weak, b 
five times stronger) to a linear bare band. 



strong electron-electron and electron-phonon interac- 
tions. Nonetheless, in the superconducting state of the 
copper-oxides at optimal and overdoped regime, one re- 
covers narrow peaks (20~30meV) of the spectral func- 
tion. The above self-energy, then, is able to describe the 
phenomenology of the mode-coupling behavior for the 
superconducting state. The difference between the self- 
energy induced for a particular mode and coupling con- 
stant in the normal state at T=100K fFig ll9() and the 
superconducting state at T=10K (FigEJ) also shows the 
extent to which one can expect a temperature-dependent 
mode coupling in the high-Tc cuprates. 

To illustrate the salient features of mode-coupling on 
the dispersion, we consider a linear bare band coupled to 
an Einstein phonon in the normal state at T=10K. The 
effect of electron-phonon interaction on the one-electron 
spectral weight A{k,uj) of a d^2_„ 2 sup erconductor has 
been simulated by Sandvik et al. l)l22f) . In Fig l20l we 
show image plots, EDCs, MDC derived dispersions, and 
the MDC extracted widths for two different coupling con- 
stants (the case of stronger coupling is a factor of five 
increase in the vertex, g^). 

There are three characteristic signatures of mode cou- 
pling behavior evident: 

1) A break up of a single dispersing peak into two 



branches (Figl20bl and bl) — a peak that decays as it 
asymptotically approaches the mode energy (I in Fig. 
I20b.2 and b2), and a hump that traces out a dispersing 
band (II in Fig|2nt2 and b2). 

2) In the image plots fFigl20bl and bl), a significant 
broadening of the spectra beyond the mode energy is 
readily apparent. This is also the origin of the broad 
hump of the dispersing band seen in the EDCs (Figl20b2 
and b2) and the step in the extracted widths (or lifetime) 
(Fig|23i4 and b4). 

3) At the mode energy itself, there is a "dip" between 
the peak and the hump in the EDCs (Fig. I^nk2 and b2) 
leading to the "peak-dip- hump" structure often discussed 
in the literature. 

From these generic features of electron-phonon cou- 
pling, one could ascertain the mode energy and coupling 
strength. Theoretically, the mode energy should be the 
energy to which the peak in the EDC curve decays. If 
there is a well-defined peak that has enough phase space 
range to decay, the last point at which it can be mea- 
sured is the best indication of the mode energy. Oth- 
erwise, estimates can be made from the EDC, MDC- 
derived dispersions, and the position of the step in the 
MDC widths. The coupling strength is indicated by the 
extent of the break up of the spectra into a peak and a 
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hump, the sharpness of the "kink" in the MDC-derived 
dispersion, and the magnitude of the step in the MDC- 
derived widths. Quantitative assessments of the couphng 
strength, however, require either a fuU model calculation 
or an extraction procedure to invert the phonon density 
of states coupled to the electronic spectra. 



3. Strong Coupling - Polaron 

When the kinetic energy of the particles is less than 
the phonon energy, the dressing of the phonon cloud 
could be large and the el-ph coupling enters into the po- 
laron regime. A single particle coupled to the phonon 
is the typical case, on which extensive theoretical stud- 
ies have been done. Let g{q) be the coupling con- 
stant of the phonon with wavenumber q to the electrons, 
and the lattice relaxation energy Elr is estimated as 
Elr =< \giq)\'^ > /^- When this Elr is smaller than 
the bandwidth, the effective reduction of the el-ph cou- 
pling due to the rapid motion of the electron, i.e., the 
motional narrowing, occurs and the weight of the one- 
phonon side-band is of the order of g{q)/W with the 
number of the phonon quanta iVph being estimated as 
Nph ~< |<?(9)P > /W^ - S{n/W)^ where S = ELR/n. 
As the el-ph coupling constant increases, the polaron 
state evolves from this weak coupling large polaron to 
the strong coupling small polaron. This behavior is non- 
perturbative in nature, and the theoretical analysis is 
rather difficult. One useful method is the adiabatic ap- 
proximation where the frequency of the phonon is set 
to be zero while E^n remains finite. In this limit, one 
can regard the phonon as a classical lattice displacement, 
whose Fourier component is denoted by Qq. Then one 
can investigate the stability of the weak coupling large 
polaron state, i.e., zero distortion state in the present 
approximation, by the perturbative way. Namely the en- 
ergy gain second order in g{q) reads: 
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(4.8) 



with the energy dispersion E{k) of the electron. Here the 
electron is at the ground state with the energy E{0) in 
the unperturbed state. Introducing the index i charac- 
terizing the range of the coupling as g{q) oc q^^, and con- 
sidering the smallest possible wavenumber g,„in oc N^^''^ 
for the linear sample size L — N^''^ in spatial dimension 
d, one can see that the index 



d-2{l+l) 



(4.9) 



separates the two different behavior of 5E. For s > 0, 5E 
for q = gniin goes to zero as A^ ^ oo, which suggests that 
the weak coupling state is always locally stable, separated 
by an energy barrier from the strong coupled small po- 
laron state. This means that a discontinuous change from 
the weak to strong coupling polaron states occurs where 



the mass becomes so heavy that the carrier is easily local- 
ized by impurities. Namely, the self-trapping transition 
occurs. For s < 0, on the other hand, the zero distor- 
tion state is always unstable for infinitesimal g{q) and 
hence the lowest energy state continues smoothly as the 
coupling increases, i.e., no self-trapping transition. The 
most relevant case of the short range el-ph coupling in 
two-dimensions, i.e., d = 2, £ = 0, corresponds to s = 0, 
and hence is the marginal class. Therefore whether the 
self-trapping transition occurs or not is determined by 
the model of interest, and is nontrivial. 

For the study of the polaron in the intermediate 
to strong coupling region, one needs to invent a reli- 
able theoretical method to calculate the energy, phonon 
cloud, effective mass, and the spectral function. Up to 
very recently, it has been missi ng b ut the diagrammatic 
quantum Monte Carlo method l)l23t) c ombined with the 
stochastic analytic continuation (124) enabled the "nu- 
merically exact" solution to this difficult problem. By 
this method, the crossover from the weak to strong cou- 
pling regions have been analyzed accurately for various 
models (,125; ,126) With this method, the polaron prob- 
lem in the t-J model has been studied, and detailed in- 
formation on the spectral function is now available which 
can be directly compared with experimental results. It 
is found that the self-trapping transition occurs in the 
two-dimensional t-J polaron model, and in comparison 
with experiment, the realistic coupling constant for the 
undoped case corresponds to the strong coupling region. 
Namely the single hole doped into the undoped cuprates 
is self-trapped. See below (IV. D) for more details of how 
the polaron model relates to such experimentally deter- 
minable quantities as the lineshape, dispersion, and the 
chemical potential shift with doping. 

Now we turn to the ARPES measurements that can 
be related to the two regimes of electron-phonon cou- 
pling. We will first review the band renormalization ef- 
fects along the (0,0)-(7r,7r) nodal direction and near the 
(7r,0) antinodal region. The weak electron-phonon cou- 
pling picture is useful in accounting for many observa- 
tions. However, there are experimental indications that 
defy the conventional electron-phonon coupling picture. 
Then we will move on to review the polaron issue which 
manifests in undoped and heavily underdoped samples. 



C. Band Renormalization and Quasiparticle Lifetime Effects 

1. El-Ph Coupling Along the (0,0)-(7r,7r) Nodal Direction 

The nodal direction denotes the (0,0)-(7r,7r) direction 
in the Brillouin zone (Fig. ^p). The rf-wave super- 
conducting gap is zero along this particular direction. 
As shown in Fig|2] and Figl22b. the energy-momentum 
dispersion curves from MDC method exhibit an abrupt 
slope change ("kink") near 70 meV. The kink is accompa- 
nied by an accelerated drop in the MDC width at a simi- 
lar energy scale (Fig. 122b). The existence of the kink has 
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FIG. 21 Ubiquitous existenc e of a kink in the nodal dispersion 
of various cuprate materials l|l23). Top panels (a, b, c) plot 
dispersions along (0, 0)-(7r,7r) direction (except for panel b 
inset, which is off this direction) as a function of the rescaled 
momentum k' for different samples and at different doping 
levels (5): (a) LSCO at 20 K, (b) Bi2212 in superconducting 
state at 20 K, and (c) Bi2201 in normal state at 30 K. Dotted 
lines are guides to the eye. Inset in (b) shows that the kinks 
in the dispersions off the (0, 0)-(7r,7r) direction sharpen upon 
moving away from the nodal direction. The black arrows indi- 
cate the position of the kink in the dispersions. Panels (d) and 
(e) show the temperature dependence of the dispersions for 
(d) optimally doped LSCO (x=0.15) and (e) optimally-doped 
Bi2212, respectively. Panel (f) shows the doping dependence 
of the effective electron-phonon coupling strength A' along the 
(0, 0)-(7r,7r) direction. Data are shown for LSCO (filled trian- 
gles), Nd-doped LSCO (1/8 doping; filled diamonds), Bi2201 
(filled squares), and Bi2212 (filled circles in the first Brillouin 
zone and unfilled circles in the second zone). The different 
shadings represent data obtained in different experimental 
runs. Shaded area is a guide to the eye. 



been well established as ubiquitous in hole-doped cuprate 
materials(127; 128; 129; 130; 131; 132; 13l): 

1. It is present in various hole-doped cuprate mate- 
rials, including BiaSraCaCusOg (Bi2212), BisSrsCuOg 
(Bi2201), (La2-a;Sra;)Cu04 (LSCO) and others. The en- 
ergy scale (in the range of 50-70 meV) at which the kink 
occurs is similar for various systems. 

2. It is present both below T^ and above Tc- 

3. It is present over an entire doping range (Fig. 122b). 
The kink effect is stronger in the underdoped region and 
gets weaker with increasing doping. 

While there is a consensus on the data, the exact mean- 
ing of the data is still under discussion. The first issue 
concerns whether the kink in the normal state is related 
to an energy scale. Valla et al. argued that the system 
is quantum critical and thus has no energy scale, even 
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FIG. 22 Doping dependence of the nodal electron dynamics 
in LSCO and universal nodal Fermi velocitv<ll32D. (a) Dis- 
persion of LSCO with various doping levels (x=0.03 to 0.30) 
measured at 20 K along the (0,0)-(7r,7r) nodal direction. The 
arrow indicates the position of kink that separates the dis- 
persion into high-energy and low-energy parts with different 
slopes, (b). Scattering rate as measured by MDC width (fuU- 
width-at-half-maximum, FWHM) of the LSCO (x = 0.063) 
sample measured at 20 K. The arrow indicates a drop at an 
energy ~70 meV. 



though a band renormalization is present in the data dSJI) . 
Since their data do not show a sudden change in the scat- 
tering rate at the corresponding energy, they attributed 
the kink in Bi2212 above Tc to the margin al Fe rmi liquid 
(MFL) behavior without an energy scale l)l30|) . Others 
believe the existence of energy scale in the normal and 
superconducting states has a common origin, i.e., cou- 
pling of quasiparticles with low-energy collective excita- 
tions (bosons) (,127: 128: 129) . The sharp kink structure 
in dispersion and concomitant existence of a drop in the 
scattering rate which is becoming increasingly clear with 
the improvement of signal to noise in the data, as ex- 
emplified in und erdop ed LSCO (x=0.063) in the normal 
state (Fig. I22b) l)l34|) . are apparently hard to reconcile 
with the MFL behavior. 

The existence of a well-defined energy scale over 
an extend ed te mperature range is best seen in Bi2201 
compound l)l35tl . As shown in Fig[531 the spectra reveal 
a "peak-dip-hump" structure up to temperatures near 
130K, almost ten times the superconducting critical tem- 
perature Tc. Such a "peak-dip- hump" structure is very 
natural in an electron-phonon coupled system, but will 
not be there if there is no energy scale in the problem as 
argued by Valla et al. lJ33) . 

A further issue concerns the origin of the bosons 
invol ved in the coupling, with a m agne tic resonance 
mode l)l29l ; Il30l) and optical phonons l)l28|) being possi- 
ble candidates considered. The phonon interpretation is 
based on the fact that the sudden band renormalization 
(or "kink") effect is seen for different cuprate materi- 
als, at different temperatures, an d at a similar energy 
scale over the entire doping range l)l28(l . For the nodal 
kink, the phonon considered in the early work was the 
half-breathing mode, which shows an anomaly in neu- 
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FIG. 23 Energy Distribution curves (EDC) in the normal 
state of underdoped Bi22 01 (T e= lOK) at several tempera- 
tures (from 20K to 200K)Jl^. A dip in the EDCs can be 
clearly observed almost for all the temperatures. The dip 
position (dotted line) is 60meV and is roughly temperature 
independent. 



tron exDeriments ljlOTt Il08(l . Unlike the phonons, which 
are similar in all cuprates, the magnetic resonance (at 
correct energy) is observed only in certain materials and 
only below Tc. The absence of the magnetic mode in 
LSCO and the appearance of magnetic mode only be- 
low Tc in some cuprate materials are not consistent with 
its being the cause of the universal presence of the kink 
effect. Whether the magnetic resonance can cause an y 
additional effect is still an active research toDic l)l36Hl37|) . 

To test the idea of electron-phonon coupling, a n iso - 
tope exchange experiment has been carried out l)l33|) . 
When exchanging ^^O and ^®0 in Bi2212, a strong iso- 
tope effect has been reported in the nodal dispersions 
(Fig l24|l . Surprisingly, however, the isotope effect mainly 
appears in the high binding energy region above the kink 
energy; at the lower binding energy near the Fermi level, 
the effect is minimal. This is quite different from the con- 
ventional electron-phonon coupling where isotope substi- 
tution will result in a small shift of phonon energy while 
keeping most of the dispersion intact. The origin of this 
behavior is still being investigated. 

It is interesting to note in Fig. 122b that the energy 
scale of the kink also serves as a dividing point where the 
high and lo w en ergy dispersions display different doping 
dependence(|l32l). The dispersion in this Figure were ob- 
tained by the MDC method. In Fig. 125b .. we reproduce 
some of these MDC-extracted dispersions, but we also 
plot the dispersion extracted using EDCs by following the 
EDC peak position. Since the first derivative of the dis- 
persion, 9E/9k, corresponds to velocity, the dispersions 
at high binding energy (-0.l~-0.25eV) and low binding 
energy (0~-0.05eV) are fitted by straight lines to quan- 
titatively extract velocities, as plotted in Fig. E5b l|l38h . 

While nodal data clearly reveal the presence of cou- 
pling to collective modes with well-defined energy scale, 
there are a couple of peculiar behaviors associated with 
the doping evolution of the nodal dispersion. One ob- 
vious anomaly is the difference of low energy velocity 
obtained from MDC and EDC methods(Fig. ^^). As 
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FIG. 24 I sotope-induced changes of the nodal 
dispersion |(133[) . The data were taken on optimally 

doped Bi2Sr2CaCu208 samples (Tc ~ 91 to 92 K) with 
different oxygen isotopes ^®0 and ^^O at T ~ 25 K along 
the nodal direction. The low energy dispersion is nearly 
isotope-independent, while the high energy dispersion is 
isotope-dependent. The effect is reversible by isotope 
re-substitution (green). Inset shows the real part of the 
electron self-energy, ReE , obtained from the dispersion by 
subtracting a line approximation for the one-electron band 
Ek, connecting two points (one at Ef and the other at a 
300-meV binding energy) of the ^*0 dispersion. 



seen from Fig. 122b and Fig. 125b . the low-energy dis- 
persion and velocity from the MDC method is insensi- 
tive to doping over the entire doping range, giving the 
so-called "universal nodal Fermi velocity" be havio r! 132). 
Similar behavior was also reported in Bi2212 l)l3C]() . How- 
ever, improved LSCO data where we can resolve a well- 
defined quasiparticle peak to extract dispersion using 
EDC method reveal a dichotomy in EDC and MDC de- 
rived disper sions , particularly for low doping (Fig. I25|) . 
like x=0.01 l)l39(l . This discrepancy between EDC and 
MDC cannot be reconciled within the conventional el- 
ph interaction picture, as simulations considering exper- 
imental resolutions show. 

In terms of conventional electron-phonon coupling, if 
one considers that the "bare band" does not change 
with doping but the electron-phonon coupling strength 
increases with decreasing doping, as it is probably the 
case for LSCO, one would expect that the low energy 
dispersion and velocity show strong doping dependence, 
while the high-energy ones converge. However, one sees 
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FIG. 25 (a). Energy- momentum dispersions for LSCO wit h 
different dopings, using both EDC and MDC methods lll38fl . 
The EDC low-energy velocity is obtained by fitting the EDC 
dispersion linearly in the intermediate energy range because 
the data points very close to Fermi level is affected by the 
Fermi cutoff while the data at higher energy have large uncer- 
tainty because the EDCs are broader. The MDC low (high) 
energy velocity yh is obtained by fitting MDC dispersion at 
binding energy 0~50meV (100 250meV) using a linear line. 
(b). Low and high-energy velocities as a function of doping 
obtained from MDC and EDC dispersions. 



that the high-energy velocity is highly doping dependent. 
Moreover, its trend is anomalous if one takes electron- 
electron interaction into account. It is known in cuprates 
that, with decreasing doping, the electron-electron inter- 
action gets stronger. According to conventional wisdom, 
this would result in a larger effective mass and smaller 
velocity. However, the doping dependence of the high- 
energy velocity is just opposite to this expectation, as 
seen from Fig. 125b . 

Therefore, these anomalies indicate a potential devi- 
ation from the standard Migdal-Eliashberg theory and 
the possibility of a complex interplay between electron- 
electron and electron-phonon interactions. As we dis- 
cuss later, this phenomenon is a hint of polaronic effect 
where the traditional analysis fails. Such a polaron effect 
gets stronger in deeply underdoped system even along the 
nodal direction. Under such a condition, one needs to use 
EDC derived dispersion when the peaks are resolved. 



2. Multiple Modes in the Electron Self-Energy 

In conventional superconductors, the successful ex- 
traction of the phonon spectral function, or the Eliash- 
berg function, a^i^(w), from electron tunneling data 
played a decisive role in cementing the consensus on th e 
phonon-mediated mechanism of superconductivity (jM^)- 
For high temperature superconductors, the extraction of 



the bosonic spectral function can provide fingerprints for 
more definitive identification of the nature of bosons in- 
volved in the coupling. 

In principle, the ability to directly measure the disper- 
sion, and therefore, the electron self-energy, would make 
ARPES the most direct way of extracting the bosonic 
spectral function. This is because, in metals, the real part 
of the electron self-energy ReE is related to the Eliash- 
berg function a^F(il;e,k) by: 



ReE(k, e; T) 
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with /(x) being the Fermi distribution function. Such 
a relation can be extended to any electron-boson cou- 
pling system and the function a^-F(w) then describes the 
underlying bosonic spectral function. We note that the 
form of Rel](k,e;r) (Ea. l4.1U|l can be derived by tak- 
ing the Kramers-Kronig transformation of ImE for the 
normal state as shown above (Eg. 14. 71) . Unfortunately, 
given that the experimental data inevitably have noise, 
the traditional least-square method to invert an integral 
problem is mathematically unstable. 

Very recently, Shi et al. have made an impor- 
tant advance in extracting the Eliashberg function 
from ARPES data by employing the maximum entropy 
method (MEM) and s uccessfully applied the method to 
Be surface states lll4l|) . The MEM approach l)l4lj) is ad- 
vantageous over the least squares method in that: (i) It 
treats the bosonic spectral function to be extracted as a 
probability function and tries to obtain the most prob- 
able one. (ii) More importantly, it is a natural way to 
incorporate the priori knowledge as a constraint into the 
fitting process. In practice, to achieve an unbiased inter- 
pretation of data, only a few basic physical constraints to 
the bosonic spectral function are imposed: (a) It is pos- 
itive, (b) It vanishes at the limit uj^O. (c) It vanishes 
above the maximum energy of the self-energy features. 
As shown in the case of Be surface state, t his m ethod is 
robust in extracting the Eliasberg function l)l4l|) . 

Initial efforts have been made to extend this approach 
to underdoped LSCO and evidence for electron coupling 
to several phonon modes has been revealed l)l42|) . As seen 
from Fig. I^Hl from both the electron self-energy(Fig. 
126b.). and the derivative of their fitted curves ((Fig. 126b.). 
one can identify two dominant features near '^ 40 meV 
and ^60 meV. In addition, two addition mod es may also 
be present near ~25meV and ^75 meV ()l42|) . The mul- 
tiple features in Fig. 126b show marked difference from 
the magnetic excitation spectra measured in L SCO which 
is mostly featureless and doping dependent (|14^. In 
comparison, they show more resemblance to the phonon 
density-of-states (DOS), measured from neutron scatter- 
ing on LSCO (Figl26b) (,144,) . in the sense of the number of 
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modes and their positions. This similarity between the 
extracted fine structure and the measured phonon fea- 
tures favors phonons as the bosons couphng to the elec- 
trons. In this case, in addition to the half-breathing mode 
at 70^80 meV tha t we previously considered strongly 
coupled to electrons l|l28l) . the present results suggest that 
several lower energy optical phonons of oxygens are also 
actively involved. Particularly we note that the mode at 
'^ BOmeV corresponds to the vibration of apical oxygens. 

We note that, in order to be able to identify fine 
structure in the electron self-energy, it is imperat ive t o 
have both high energy resolution and high statistics l)l45fl . 
These requirements have made the experiment highly 
challenging because of the necessity to compromise be- 
tween two conflicting requirements for the synchrotron 
light source: high energy resolution and high photon flux. 
Further improvements in photoemission experiments will 
likely enable a detailed understanding of the boson modes 
coupled to electrons, and provide critical information for 
the pairing mechanism. 

One would like to extend this method to the super- 
conducting state, in momentum around the BZ, and to 
higher temperatures. The superconducting state could, 
in principle, be achieved by using the BCS dispersion 
of the quasiparticles rather than the normal state dis- 
persion and is currently under study. However, consid- 
ering the anisotropy of the el-ph coupling detailed be- 
low, the anisotropy of the underlying band structure, 
and the d-wave superconducting gap, extending the pro- 
cedure in momentum may be somewhat more difficult. 
The a^F^uj, e, k) used for the above form of the real part 
of the self-energy is assumed to be only weakly dependent 
on the initial energy e and momentum k of the electron. 
But again, one in principle could begin to consider a dif- 
ferent form of the calculated ReS and then apply the 
MEM method with it instead. Extending the method 
to higher temperatures, for example ^ lOOK for normal 
state Bi2212 data, may be, however, a limitation that 
cannot be overcome. The method's strength is in resolv- 
ing fine structures due to the phonon density of states. 
Those fine structures occur predominantly at lower tem- 
peratures. At higher temperatures of ^ lOOK, the imag- 
inary and real parts of the self energy get broadened on 
the order of the phonon energy itself. In that case, two 
or more neighboring phonons would contribute to the 
electronic renormalization at a given energy, both broad- 
ening the fine structures in the data and weakening the 
resolving power of the method itself. So, while the MEM 
method can directly extract fine features from ARPES 
data in agreement with neutron scattering without im- 
plicitly assuming a phonon model, it does not have the 
freedom to incorporate the temperature and momentum 
dependence needed to describe the ARPES data in both 
superconducting and normal states, near the vHS and 
near the node. Both modelling of the data and direct 
extraction, then, are needed, to gain a complete picture 
of the mode-coupling features in the data. 
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FIG. 26 M ultiple modes coupling in electron self-energy in 
LSCO (|l42t) . (a). The effective real part of the electron self- 
energy for LSCO x=0.03 (al), 0.063 (a2), 0.07 (a3)and ~0.06 
(a4) samples. Data (al-a3) were taken using Scienta 2002 
analyzer while data (a4) were taken using Scienta R4000 ana- 
lyzer. Data (al-a3) were taking using lOeV pass energy at 
an energy resolution of ~18meV. Data (a4) were taken a 
X ~0.06 sample using 5eV pass energy with an energy res- 
olution of ~12meV. For clarity, the error bar is only shown 
for data (a4) which becomes larger with increasing binding 
energy. The arrows in the figure mark possible fine structures 
in the self-energy. The data are fitted using the maximum 
entropy method (solid red lines). The values of (ai, 02) (the 
unit of OLi and a.2 are eV-A and eV-A^, respectively) for bare 
band are (-4.25,0) for (al), (-4.25, 13) for (a2), (-3.7,7) for 
(a3) and (-4.3, 0) for (a4). (b). The second-order derivative 
of the calculated ReE. The ruggedness in the curves is due 
to limited discrete data points. The four shaded areas corre- 
spond to energies of (23~29), (40~46), (58~63) and (75~85) 
meV where the fine features fall in. ((c) The phonon den- 
sity of state F{uj) for LSCO x = (red) and x = 0.08 (blue) 
measured from neutron scattering |l44tl . 



3. El-Ph Coupling Near the (tt.O) Antinodal Region 

The antinodal region refers to the (7r,0) region in the 
Brillouin zone where the d-wave superconducting gap has 
a maximum (Fig. Eb)- Recently, a low-energy kink 
was also identified near the (7r,0) antinodal region in 
Bi2212(^ J^a Jia 147) . This observation was made 
possible thanks to the clear resolution of the bi-laycr 
sphtting(^2g; |2^ Isol) . As there are two adjacent Cu02 
planes in a unit cell of Bi2212, these give rise to two Fermi 
surface sheets from the higher-binding-energy bonding 
band (B) (thick red curves in Fig. 127b ) and the lower- 
binding- energy antibonding band (A) (thick black curves 
in Fig. EZt). 

Consider a cut along (7r,7r)-(-7r,7r) across (7r,0) in 
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FIG. 27 Antinodal kink near (7r,0) in a heavily overdoped Bi2212 sample (Tc~58 K). (a) Normal-state data (T=85 K) near 
the antinodal region, (b) Superconducting-state data from the same sample at 10 K, showing the emergence of a dispersion 
kink in the bilayer split-B band. The B band dispersions (red curves) were determined by fitting MDC peak positions. The 
black dots represent A band EDC peak positions, (c) Brillouin zone with bonding band B (thick red) and antibonding band A 
(thick black) Fermi surfaces, as well as momentum-cut locations for panels (a) and (b) (blue bars). The two sets of thin curves 
are replicas of Fermi surface originating from the superstructure in Bi2212. 



61221^ both above T^ (Fig. EIJi) and below T^ (Fig. 
I27b)l;!5#. Superimposed are the dispersion of the bond- 
ing band determined from the MDC (red lines) and 
antibonding band from the EDC (black dots). When 
the bandwidth is narrow, the applicability of the MDC 
method in obtaining dispersion becomes questionable so 
one has to resort to the traditional EDC method. In the 
normal state, the bonding-band dispersion (Fig. 177k ) 
is nearly linear and featureless in the energy range of 
interest. Upon cooling to 10 K (Fig. I27b'l. the disper- 
sion, as well as the near-ii^i;' spectral weight, is radically 
changed. In addition to the opening of a superconduct- 
ing gap, there is a clear kink in the dispersion around 40 
meV. 

Gromko et al. l|54J) reported that the antinodal kink 
effect appears only in the superconducting state and 
gets stronger with decreasing temperature. Their 
momentum-dependence measurements show that the 
kink effect is strong near (7r,0) and weakens dramati- 
cally when the momentum moves away from the (7r,0) 
point. Excluding the possibility that this is a by-product 
of a superconducting-gap opening, they attributed the 
antinodal kink to the coupling of electrons to a bosonic 
excitation, such as a phonon or a collective magnetic ex- 
citation. The prime candidate they considered is the 
magnetic-resonance mode observed in inelastic neutron 
scattering experiments. 

The temperature and momentum dependence identi- 
fied for a range of doping levels has also led others t o 
attribute the effect to the magnetic resonance l|l29t Il46^ . 
However, there are some inconsistencies with this inter- 
pretation: (1) the magnetic resonance has not yet been 



observed by neutron scattering in such a heavily doped 
cuprate and (2) the magnetic resonance has little spec- 
tral weight, and may be too weak to cause the effect 
seen by ARPES. Furthermore, the electron-phonon cou- 
pling in the early tunneling spectra, such as Pb, appeared 
prominently only in the superconducting state. The lin- 
ear MDC-derived dispersion in the normal state of Bi2212 
at (pi,0) that Gromko et. al. reports! 54) is not conclu- 
sive enough proof that the same mode does not couple 
to the electrons in the normal state. On the other hand, 
the clear determination of mode-coupling by Gromko et. 
al. in the anti-nodal region, where the gap is maximum, 
without the complication of bilayer splitting or super- 
structure, suggests that the renormalization effects seen 
by ARPES in the cuprates may indeed by related to the 
microscopic mechanism of superconductivity. 

Cuk et. al. lll47l) and Devereaux et. al. ()l48j) have re- 
cently proposed a new interpretation of the renormaliza- 
tion effects seen in Bi2212. Specifically, the key observa- 
tion that prompted them to rule out the magnetic reso- 
nance interpretation is the unraveling of the existence of 
the antinodal kink even in the normal state. This obser- 
vation was made possible by utilizing the EDC method 
because the MDC method is not appropriate when the 
assumed linear approximation of the bare band fails near 
(7r,0) where the band bottom is close to Fermi level Ep. 
Figs. I28bl. [^Sbl. and 128b show dispersions in the nor- 
mal state of an optimally-doped sample which consis- 
tently reveal a 40 meV energy scale that has eluded 
detection previously. Upon entering the superconduct- 
ing state, the energy scale shifts to 70meV consistent 
with a gap opening of 35'--^40 meV. This coupling is also 
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found to be more extended in a Brillouin zone than pre- 
viously report ed ljS^ . In Fig. 1201 we show data from the 
optimally-doped Bi2212 sample for a large portion of the 
BZ in the superconducting state(l47j. The renormaliza- 
tion occurs at 70 mcV throughout the BZ and increases 
in strength from the nodal to anti-nodal points. Similar 
behaviors are also noted by others l)l29fl (Fig. I31|) . The 
increase in coupling strength can be seen in the follow- 
ing ways: Near (7r,0), the band breaks up into two bands 
(peak and hump) as seen in Fig. I30b2 and a3. For cuts 
taken in the (0,0) - (tt, it) direction, the band disper- 
sion is steeper and the effects of mode-coupling, though 
significant, are less pronounced. 

It is quite natural that phonon modes of different origin 
and energy preferentially couple to electrons in certain 
k-space regions. While the detection of multiple modes 
in the normal state of LSCO((142) suggests that sev- 
eral phonons may be involved, only one has the correct 
energy and momentum dependence to understand the 
prominent signature se en in the superconducting state. 
This new interpretation l|147|) attributes the renormaliza- 
tion seen in the superconducting state to the "bond- 
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FIG. 28 Antinodal kink near (7r,0) in the normal state 
(al,bl,c) and superconducting state (a2,b2) in an optimally- 
doped Bi2212|J5J). The dispersions in al, bl, and c were 
derived by the EDC method; the position of the momentum 
cuts is labelled in the insets. The red dots are the data; the 
fit to the curve (black dashed line) below the 40-meV line 
is a guide to the eye. The dispersions at the same location 
in the superconducting state (10 K) are shown in (a2) and 
(b2), which were derived by the MDC method (blue circles). 
Also plotted in (a2) and (b2) are the peak (blue squares, I) 
and hump positions (blue squares, II) of the EDCs for com- 
parison. The inset of (a2) shows the expected behavior of 
a Bogoliubov-type gap opening. The s-like shape below the 
gap energy is an artifact of the way the MDC method han- 
dles the backbend of the Bogoliubov quasiparticle. (d). Kink 
positions as a function of momentum cuts in the antinodal 
region. 



buckling" Big phonon mode involving the out-of-plane, 
out-of-phase motion of the in- plane oxygens. The bond- 
buckling phonon is observed at 35 meV in the Bi^ po- 
larization of Raman scattering on an optimally doped 
sample, the same channel in which the ~ 35-40 meV 
d-wave superconducting gap shows up 1)99 : Il49t Il50t) . 
Applying simple symmetry considerations and kinematic 
constraints, it is found that this Big buckling mode in- 
volves small momentum transfers an d co uples strongly 
to electronic states near the antinode(ll43)- In contrast, 
the in-plane Cu-0 breathing modes involve large momen- 
tum transfers and couple strongly to nodal electronic 
states. Band renormalization effects are also found to 
be strongest in the superconducting state near the antin- 
ode, in full agreement with angle-resolved photoemission 
spectroscopy data (Fig. I29|l . The dramatic temperature 
dependence stems from a substantial change in the elec- 
tronic occupation distribution and the opening of the su- 
perconducting gap( 147; 148). It is important to note that 
the electron-phonon coupling, especially the one with Big 
phonon, explains the temperature and momentum depen- 
dence of the self-energy effects that were taken as key ev- 
idence to support the magnetic resonance interpretation 
of the data. Compounded with the findings that cannot 
be explained by the magnetic resonance as discussed ear- 
lier, this development makes the phonon interpretation 
of the kink effect self-contained. 



4. Anisotropic El-Ph Coupling 

The full Migdal-Eliashberg-based calculation consists 
of a tight-binding band structure and el-ph coupling to 
the breathing mode as well as the Big bo nd-bu ckling 
mode and is based on an earlier calculation {15¥). The 
electron-phonon coupling vertex g{k,q), where k repre- 
sents the initial momentum of the electron and q the mo- 
mentum of the phonon is determined on the basis of the 
oxygen displacements for each mode in the presence of 
the underlying band-structure. In the case of the breath- 
ing mode, the in-plane displacements of the oxygen mod- 
ulate the Cu02 nearest neighbor hopping integral as well 
as the site energies. In the case of the bond-buckling 
mode, one must suppose that the mirror plane symmetry 
across the Cu02 plane is broken in order for electrons to 
couple linearly to phonons. The mirror plane symmetry 
can be broken by the presence of a crystal field perpen- 
dicular to the plane, tilting of the Cu-0 octahedral, static 
in-plane buckling, or may be dynamically generated. 

The gBi (k, q) form factor leads to preferential q ~ 2k f 
scattering between the parallel pieces of Fermi surface 
in the anti-nodal region, as shown in Fig. 1321 depicting 
g(fc, k') for the buckling mode (where k' — k ~ q) for an 
electron initially at the anti-node [kAN] upper-left) and 
for an electron initially at the node (fc^r; bottom- left) . 
This coupling anisotropy partially accounts for the strong 
manifestation of electron-phonon coupling in the anti- 
nodal region where one sees a break up into two bands. 
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FIG. 29 Comparison between the calculated and measured spectral function in Bi2212 including electron-phonon coupling for 
three different momentum cuts (a, b, c) through the Brillouin zone . (al,bl,cl) and (a2,b2,c2) show the calculated spectral 
functions in the normal and superconducting states, respectively i 14^ . The measured spectral functions are shown in (a3,b3,c3) 
for the normal state and in (a4,b4,c4) for the superconducting state. The corresponding momentum cuts a, b, and c are shown 
in the rightmost panel. The red markers in the superconducting state indicate 70 meV. The simulation includes Big oxygen 
buckling mode and half-breathing mode. 
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FIG. 30 Anisotropic electron-boson coupling in Bi2212((l43). Image plots in (al-a6) and (bl-b6) are cuts taken parallel to (0, 
7r)-(7r, n) and (0, 0)-(7r, n) respectively at the locations indicated in the Brillouin zone ((a) and (b)) at 15 K for an optimally 
doped sample (Tc = 94 K). 



The breathing mode, in contrast, modulates the hopping 
integral and has a form factor, gbr{k,q), that leads to 
preferential scattering for large q and couples opposing 
nodal states. This coupling anisotropy then accounts for 
the 70 meV energy scale seen most prominently in a nar- 
row /c-space region near the nodal direction in the normal 
state of LSCO. Fig. |221also shows that the magnitude 
of the electron-phonon vertex is largest for an electron 
initially sitting at the node, k^, that scatters to the op- 
posing nodal state. F or mo re details on this calculation, 
see Devereaux et. al. lll48l). 



The anisotropy of the mode-coupling in both the su- 
perconducting state data and the calculation is pecu- 
liar to the cuprates. Such a strong anisotropy in the 
electron-phonon coupling is not traditionally expected. 



In cuprates, the sources of the anisotropy are: 1) an 
electron-phonon vertex for the Big bond-buckling mode 
and the breathing mode that depends both on the elec- 
tron momentum k as well as the phonon momentum q. 
This comes from a preferential scattering in the Brillouin 
zone, in which nodal states couple to other nodal states 
and anti-nodal states to other anti-nodal states. 2) a 
strongly anisotropic electronic band structure character- 
ized by a van Hove Singularity (vHS) at (7r,0). In the 
anti-nodal region and along the (tt, 0) — (tt, tt) direction in 
which 2fci? scattering is preferred, the bands are narrow, 
giving rise to a larger electronic density of states near the 
phonon energy and therefore a stronger manifestation of 
electron-phonon coupling. 3) a d-wave superconducting 
gap and 4) a collusion of energy scales in the anti-nodal 
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FIG. 32 Plots of the electron-phonon coupling | g{]i,q) p for 
initial k and scattered k' = k — q states on the Fermi surface 
for the buckling mode (left panels) and breathing mode (right 
panels) for initial fermion k at an anti-nodal (top panels) 
and nodal (bottom panels) point on the Fermi surface, as 
indicated by the arrows. The red/blue color indicates the 
maximum /mini mum of the el-ph coupling vertex in the BZ 
for each phonon( 148t) . 



FIG. 31 Momentum de pend ence of photoemission data in 
optimally-doped Bi2212 il29l) . Left panels: Photoemission 
data in the normal state (T=140 K) along selected cuts par- 
allel to M(7r,0)-Y(7r,7r). EDC peak positions are indicated by 
crosses. Middle panels: Photoemission data in the supercon- 
ducting state (T=40 K) at the same cuts as for left panels. 
Crosses indicate positions of broad high energy peaks, dots 
sharp low energy peaks. Right panels: EDCs at locations 
marked by the vertical lines in the middle panels. 



region that resonate to enhance the above effects — the 
vHS at r^ SSmeV in the tight-binding model that best 
fits the data, the maximum d-wave gap at ~ 35meV, 
and the bond-buckling phonon energy at ~ 35meV. All 
these three factors collide to give the anisotropy of the 
mode-coupling behavior in the superconducting and nor- 
mal states. For a detailed look at how each plays a role 
in th e agreement with the data, please see Cuk et. al. 
l)l47fl . The coincidence of energy scales, along with the 
dominance of the renormalization near the anti-node, in- 
dicates the potential importance of the Big phonon to the 
pairing mechanism, wh ich is consisten t with some theory 
on the Big phonon l(89t Il52t Il53t Il54|) but remains to be 
investigated. 

The cuprates provide an excellent platform on which 
to study anisotropic electron-phonon interaction. In one 
material, such as optimally doped Bi2212, the effective 
couphng can span A of order ~1 at the node to 3 at the 
anti-node (Fig lH5|l l|l47l: Hii). In addition to the large 
variation of coupling strength, there is a strong varia- 




FIG. 33 Plots of the electron-phonon coupling Ak in the first 
quadrant of the BZ for the buckling mode (right top panel) 
and breathing mode (right bottom panel). The color scale is 
shown on the right for each phonon. The left panel shows 
energy contours for the band structure used jlMSJ) . 



tion in the kinematic considerations for electron-phonon 
coupling. In the nodal direction, the band bottom is 
far from the relevant phonon energy scales. However, at 
the anti-node, the relevant phonon frequencies approach 
the bandwidth. Indeed the approximation of Migdal, in 
which higher order vertex corrections to the el-ph cou- 
pling are neglected due to the smallness of (A * fi i^/Ep), 
may be breaking down in the anti-nodal region. Non- 
adiabatic effect beyond the Migdal approximation have 
been considered and are under continuing study (^155(1 . 
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FIG. 34 (a). Photoemissi on spectr a along the (0,0 ) and (tt.tt ) 
direction in Sr7Cu07Cl9 l|l5fil: \lE^ Ca oCuOoClo llTTl llfil). 
La2Cu04(|l6i[l63) and Nd^CuO^Jiel. 



D. Polaronic Behavior 

1. Polaronic Behavior in Parent Compounds 

The parent compounds of the cuprate superconduc- 
tors, being antiferromagnetic Mott insulators, provide 
an ideal testing ground for investigating the dynamics 
of one hole in an antiferromagnetic background. In- 
deed, many theories have been formed and tested by 
ARPES on a number of compounds, among them are 
Sr^ Cu09Cl9a56; 157: 158: 159. 160) , CaoCuOoCL([Til: 

166l : ll67J) . However, several aspects of the data can only 



be explained by invoking polaron physics, as we will now 

discuss. 

The ARPES measurements on SCOCljlli [HI) and 
CCOC l)ll3Hl62|) give essentially similar results. As seen 
in Fig. 134b , and 134b . along the (0,0)-(7r,7r) direction, 
the lowest energy feature disperses toward lower bind- 
ing energy with increasing momentum, reaches its lowest 
binding energy position near (7r/2,7r/2) where it becomes 
sharpest in its lineshape, and then suddenly loses inten- 
sity after passing (7r/2,7r/2) and bends back to high bind- 
ing energy. This behavior can be more clearly seen in the 
image plot of Fig. I35b ( 162) where the "band" breaks into 
two branches. The lowest binding energy feature shows a 
dispersion of ^^0.35 eV while an additional band at high 
binding energy (Fig. 135b ) is very clo se to the unrenormal- 
ized band predicted by band theory! 162). The dispersion 
of low binding-energy band along the (0,0)-(7r,7r) direc- 
tion and other high symmetry directions are shown in 
Fig. 135b by keeping track of the EDC peak position(51j. 
The total dispersion of the peak is ~0.35eV. This is in 
contrast to the predictions of one-electron band calcula- 
tions which gives an occupied band width of ^ 1.5 eV and 
total bandwidth of ~3.5 eV(168). Nevertheless, it is con- 
sistent with the calculations from the t- J mod el wh ere the 
predicted occupied bandwidth is ^2.2J l|5(]t Il69l) . This 
indicates that the dynamics of one-hole in an antiferro- 
magnetic background is renormalized from scale t to scale 
J. 
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FIG. 35 (a). Intensity plot of ARPES data as functions of 
the binding energy and mome ntum for Ca2Cu02Cl2 along 
the r(0,0)-(7r,7r) direct ion |(l^2)- The data was symmetrized 
around the F point. Also shown on the plot are the disper- 
sions obtained by following the peak positions of the MDCs 
(solid line) and the EDCs (circle and triangles). The results 
are compared with the shifted dispersion from the LDA cal- 
culation (dashed line), (b) Energy dispersion of quasiparti- 
cle for insulating Sr2Cu02Cl2 measured from the top of the 
band. Experimental data a re ta ken from (11561) (open circles), 
lll57t) ('open triangles) and il58D fopen squares). Solid circles: 
the results of the self-consistent Born approximation (SCBA) 
for the t-t'-t"-J model with t = 0.35 eV, t = -0.12 eV, 
t = 0.08 eV and J — 0.14 eV. The solid lines are obtained 
by fitting the SCBA data to a dispersion relation given by 



Eo(k) + Ei(k), being t.j-j- = -0.038 eV and t, 



eff 



0.022 eV. 



The dashed line along the (tt, )C(0,7r) direction represents 
the spinon dispersion from il72l) . 



While the t-J model and experiments show agree- 
ment along the (0,0)-(7r,7r) direction, there are discrep- 
ancies along other directions, such as the (0,0)-(7r,0) and 
(0,7r)-(7r/2,7r/2)-(7r,0) directions (156). The later inten- 
sive theoretical effort resolved this issue by incorporating 
the hoppi ng t o the second (t ) and third (t") nearest- 
neighbors! ITQ). More precise calculations of the disper- 
sion in the t— t — t —J model are performe d by using a 
self-consistent Born approximation fSCBA)l )l7ll) . These 
calculations show a satisfactory agreement with experi- 
mentally derived dispersion, as shown in Fig. I35b f5l|h. 

However, there remain a few prominent puzzles re- 
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FIG. 36 (a) Photoemision spectrum of Ca2Cu02Cl2 at 
k=(7r,7r) with fits to a Lorentzian spectral function (dashed) 
and Gaussian (red or gray)(jill^. A and B denote the 
peak maximum and the onset of spectral weight, respec- 
tively. Comparison with Fermi-liquid system Sr2Ru04 is 
shown (t hin b lack). Upper inset shows photoemission spectra 
from H2l|l73). (c) Dispersion of A and B along (0,0)-(7r,7r), 
along with experimental values for the chemical potential fj, 
(lines). 



lated to the interpretation of the photoemission data 
in undoped parent compounds (113). The first promi- 
nent issue is the linewidth of the peak near (7r/2,7r/2). 
As highhghted in Fig. 136b . the width of the sharpest 
peak near (7r/2,7r/2) is ~300 meV which is compara- 
ble w ith the entire occupied bandwidth 2J w350 meV 
||ll3). This is much broader than that from t-J model 
calculations and too broad to be considered as a coher- 
ent quasiparticle peak for which the quasiparticle peak 
is basically resolution limited, as exemplified by the data 
on Sr2Ru04 in Fig. 136b . An early attempt interpreted 
this anomalously large linewidth to additional interaction 
with a non— magnetic boson bath of excitations, such as 
phonons <)l59(l . But this interpretation meets with dif- 
ficulty in explaining little renormalization in the disper- 
sion from this "extra interaction" because dispersion and 
linewidth are closel y rel ated. A diagrammatic quantum 
Monte Carlo studv ljlTST) showed that this problem can 
be resolved by considering the polaron effect in the t-J 
model. Namely the dispersion for the center of mass of 
the spectral function obeys that of the pure t-J model, 
while the lineshape is strongly modified. The details of 
this will be given below. 

Another unresolved issue is the chemical potential /i. 
For an insulator, fi is not well defined, and may be pinned 
by surface defects or impurities and will vary between 
different samples. If one considers that the peak A in 
Fig l36b represents a quasiparticle peak, one would expect 
the chemical potential to vary anywhere above the top 
of this valence band. However, the experimental chemi- 
cal potential clearly sets a lower boun d that is ~0.45ey 
apart from the peak A (Fig. OBbl lfTT^ . Shen et al. l(TT^ 
invented a new method to determine the chemical poten- 
tial using both the energy of the non-hybridized oxygen 
orbital and the detailed line-shape of Na-CCOC. 



The resolution of these discrepancies between experi- 
ment and expectation leads to identifying polaron physics 
as responsible for the bulk of the lineshape in underdoped 
cuprates. In fact, the photoemission spectra in the under- 
doped cuprates resemble the Frank-Condon effect seen in 
photoemission spectra of molecules such as H2 (173) (inset 
of Fig. Onil, where only the "0-0" peak (filled black) 
represents the H final state with no excited vibrations 
and comprises only ^--^10% of the total intensity. In the 
solid state, this "0-0" would correspond to the quasipar- 
ticle or the coherent part of the spectral function, Acoh, 
whereas the excited states comprise the incoherent part, 
Ainc- This behavior is reminiscent of polarons, and such 
models have been invo ked in systems where strong cou- 
plings are present (174) . In this picture, in the undoped 
compound, the true QP (B) is hidden within the tail of 
spectral intensity, with a quasiparticle residue Z vanish- 
ingly small, while feature A is simply incoherent weight 
associated with shake-off excitations. 

From the viewpoint of polaron physics, the cuprates of- 
fer a unique and first opportunity to compare experimen- 
tal spectra with theory in detail. The single hole inter- 
acts both with magnons and phonons. The hole-magnon 
interaction has been successfully analy zed in terms of 
the self-consistent Born approximaiton l) 1 7 l(l . The suc- 
cess of the Born approximation results from a "satura- 
tion" effect; namely the single spin 1/2 can flip only once, 
and hence magnon clouds do not become large enough 
to induce the self-trapping transition to the small po- 
laron. On the other hand, phonon clouds can be larger 
and larger as the coupling constant g increases and can 
lead to a self-trapping transition. The t-J model coupled 
to phonons in the polaronic regime has illuminated one- 
hole dynamics in the parent compound in the following 
wav l|175(l . (1). With increasing electron-phonon cou- 
pling strength, the spectral function experiences a tran- 
sition from weak-coupling, to intermediate coupling, and 
to strong-coupling regimes. (2). In the strong coupling 
regime, the spectral function consists a ground state res- 
onance (as indicated by vertical arrows) with vanishing 
intensity and a broad peak denoted as "coherent C" , as 
shown in Figl37b. (3). The broad peak C shows strong 
momentum dependence while the lowest state is disper- 
sionless. These results are in good correspondence to 
experimental observations. The most surprising result is 
that the broad resonance has the momentum dependence 
of the t-J model without coupling to phonons (shown in 
Fig. 137b). In the Franck-Condon effect for molecules a 
similar result occurs. The center of the shake-off band 
corresponds to the hole motion in the background of the 
frozen lattice configuration, i.e. the dispersion of the 
hole remains that of the non- interacting limit, while the 
line-width broadens. A more elaborated analytic treat- 
ment of the t -J po laron model in the Franck-Condon 
approximation l)l76() successfully reproduced this Monte 
Carlo results. The calculated spectral function line-shape 
most consistent with experiment has a A = 0.9— 1.3, well 
within the strong-coupling, small-polaron regime. Re- 
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FIG. 37 (a). Calculated hole spectral func tion in ground state 
at J/t — 0.3 for different moment a ([l2a). (al) Full energy 
range for k = (7r/2,7r/2). (a2) — (a4) Low energy part for dif- 
ferent momenta. Slanted arrows show broad peaks which can 
be interpreted in ARPES spectra as "coherent" (C) and in- 
coherent (I) part. Vertical arrows indicate position of ground 
state resonance which is not seen in the vertical scale of the 
figure, (b). Dispersion of resonance energies at J/t — 0.3. 
Broad resonance (filled circles) and lowest polaron resonance 
(filled squares) at <? = 0.231125; third broad resonance (open 
circles) and lowest polaron resonance (open squares) at g=0.2. 
The solid curves are dispersions of a hole in the pure t- J model 
at J/t = 0.3. 
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FIG. 38 Polarons in La^CuOd llieTD . Calculated ARPES 
spectra for the undoped La2Cu04 system at T=0 for dif- 
ferent k normalized to the height of the phonon side band. 
The lower abscissa shows binding energies (BE) and the up- 
per abscissa the energies of the final states corresponding to 
the spectral features. The inset shows the dependence of the 
width of the phonon side band on its binding energy. The 
width of the (0, 0) spectrum is poorly defined and not shown. 



cent realistic shell model calculation|167D also concluded 
A = 1.2 for La2Cu04. 

In La2Cu04, a broad feature near -0.5eV (Fig. 
was identified as the lower Hubbard band ljlBfil : ll( 
The electron-phonon coupling strength, calculated using 
a shell model, puts La2Cu04 in the polaron regime, sim- 
ilar to Ca2Cu02Cl2. In this picture, the -0.5eV feature 
corresponds to the phonon side-band while the real quasi- 
particle residue is very weak. As shown in Fig. 1381 the 
calculated spec tral function agrees well with the mea- 
sured datalllf 



2. Doping Dependence: From Z~0 Polaron to Finite Z 
Quasiparticles 

We next turn to the question of how the small po- 
laron state evolves as a function of doping, connecting 
to the Migdal-Eliashberg regime discussed in section C. 
There are two possible ways to dope t he M ott insulator, 
schematically shown in Fig. 03; (42 Il77|) : (1). Upon 
doping, the chemical potential shifts to the top of the va- 
lence band for hole doping (Fig. llUb3) or to the bottom 
of the conduction band for electron doping (Fig. [Tnb4). 
(2). The chemical potential is pinned inside the charge 
transfer gap. Upon doping, new states will form inside 
the gap (Fig. [TO|:5). 

Recent ARPES measurements on lightly-doped 
(La2_a;Sr3,)Cu04 compounds provide a good window 
to look into this issue. As shown in Fig. 139b and 1391; . 
for undoped La2Cu04, the main feature is the broad 



peak near -0.5 eV which exhibits weak dispersion(J16fi 
There is also little spectral weight present near the 
Fermi level. However, upon only a doping of x = 0.03, 
the electronic structure undergoes a dramatic change. A 
new dispersive band near the Fermi level develops along 
the (0,0)-(7r,7r) nodal direction(Fig. 139b . right panel), 
while along the (0,0)-(7r,0) direction a saddle band 
residing -0.2eV below the fermi level develops. Even 
for more lightly-doped samples, such as cc = 0.01, new 
states near the Fermi level are created l|l39j) . Note that, 
for these lightly-doped samples, the original -0.5 eV 
remains, although with weakened spectral weight (Fig. 
I39H). So, the -0.5eV peak and the new dispersive band 
coexist at doping levels close to the parent compounds. 

The systematic evolution of the photoemission spec- 
tra near the nodal and antinodal regions with doping 
in LSCO is shown in Fig. EOK and b l)l66l) . The nodal 
quasiparticle weight, Zqp, integrated over a small en- 
ergy window near the Fermi level, is shown in Fig. 140b . 
In the underdoped region, it increases with increasing 
doping nearly linearly, and no abrupt change occurs 
near the nonsuperconductor-superconductor transition 
at X -0.05. 

(Ca2_i,Na3;)Cu02Cl2 (Na-CCOC) is another ideal 
system to address the doping evolution of the electronic 
structure. The precise measurement of the chemical po- 
tential (Fig. 141b . )■ in conjunction with the identifica- 
tion of polaron physics in the under-doped compounds, 
provides a globally con siste nt picture of the doping evo- 
lution of the cuDrates l)ll3|l . Instead of measuring the 
chemical poten tial w ith deep core level spectroscopy (the 
usual method) fl79|), one utilizes orbitals in the valence 
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FIG. 39 Creation of nodal quasiparticles in lightly-doped 
LSCO tiee') . ARPES spectra for LSCO with x=0 and j:=0.03. 
Panels a and b are EDO's along the nodal direction (0,0)- 
(tt, tt) in the second Brillouin zone (BZ). The spectra for x = 
0.03 are plotted on an expanded scale in panel c. Panels d 
and e represent energy dispersions deduced from the second 
derivative of the EDO's. 
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band at lower energies(Fig. 141b . and b). The measured 
chemical shift, A/x, exhibits a strong doping dependence, 
9/z/9x=-1.8±0.5 eV/hole, comparable to the band struc- 
ture estimation (~-1.3 eV/hole) (Fig. 1411). 

Figs. I42r a-d) show the doping evolution of the near- 
Ep EDCs plotted relative to /io of the undoped sample 
(determined in Fig. I41h). With doping, feature A evolves 
smoothly into a broad, high energy hump with a back- 
folded dispersion similar to the parent insulator (sym- 
bols), while B shifts downward relative to its position in 
the un-doped compound. Spectral weight increases with 
doping at B, and a well-defined peak emerges for the 
X =0.10 and 0.12 samples, resulting in a coherent, low- 
energy band. The dispersion of the high-energy hump 
(A) , tracked using the local maxima or second derivative 
of the EDCs, shows little change as a function of dop- 
ing (Fig. 141b). The lowest energy excitations (feature 
B, -0.05eV<E<Ei?), tracked using MDC analysis, evolve 
with doping in such a way that the quasiparticle disper- 
sion (vp) and Fermi wave vectors (kp) virtually collapse 
onto a single straight line. 



3. Doping Evolution of Fermi Surface: Nodal-Antinodal 
Dichotomy 

So far, we have discussed the doping evolution along 
the nodal direction, and seen that a sharply defined 
quasiparticle peak develops out of the small weight near 
the chemical potential in the undoped samples. We now 
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FIG. 40 ARPES spectra at k = kp in the nodal direction in 
the second BZ (a) and those at (7r,0)(b) for various doping 
levelsC166). (c). Doping dependence of the nodal QP spectral 
weight, Znqp, and the spectral weight inte grated at Ef over 
the entire second Brillouin zone, npF.R lll66tl . They show sim- 
ilar doping dependence to t he h ole concentration evaluated 
from Hall coefficient (nnaii) lIlTST) . 



FIG. 41 Chemical potential shift in Na-OCOOJTT^. (a) Va- 
lence band spectra for x=0, 0.05, 0.10, and 0.12 compositions 
at k=(0,0) and (tTjTt). O 2pz and O 2ppi states are marked by 
triangles and circles, respectively, (b) Shifts of the O 2pz and 
O 2ppi peaks shown on an expanded scale, (c) Doping de- 
pendence of chemical potential A^ determined from (b). (d) 
Valence band at k=(7r/2,7r/2), showing the lower Hubbard 
band (A) on an expanded scale. 
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FIG. 42 (a-d). EDC spectra of Na-CCOC x=0 (a), 0.05(b), 
0.10(c), and 0.12(d) from (0.27r,0.27r) to (0.727r,0.727r) with 
hump positions mark ed by open symbols and the EDC at 
kp shown in bold ljll^T) . Data are plotted on a relative energy 
scale referenced to the shift in fi shown in Fig. 141b . (e). Sum- 
mary of hump (symbols)from Fig. (a-d) and MDC dispersions 
(lines) . 



discuss the doping evolution in other directions of the 
Brillouin zone. Surprisingly, one finds that the coher- 
ent peak near the Fermi level in the lightly doped sam- 
ples is confined to the nodal region, and quickly disap- 
pears with momentum around the Brillouin zone. The 
spectral weight near the Fermi level, confined to the 
(7r/2,7r/2) nodal region, forms a so-called "Fermi arc". 
This dichotomy betw een n odal and antinodal excitations 
is shown in Fig. 1^1341) . For the x = 0.063 sample, 
which is close to the nonsuperconductor-superconductor 
transition and therefore heavily underdoped, the spec- 
tral weight near Fermi level is mainly concentrated near 
the nodal region (Fig. 143b). The coherent peaks in the 
EDCs (Fig. I43bl) near the nodal region disappear as 
one approaches the anti-nodal region, where the EDCs 
exhibit a step rather than a peak. The LSCO x = 0.09 
sample exhibits similar behavior(Fig. I43b2). In con- 
trast, for overdoped LSCO such as x = 0.22 (Fig. IIS13). 
sharp peaks are observable along the entire Fermi surface. 
These observations indicate that the electrons near the 
antinodal region experience additional scattering. There- 
fore, as shown in Fig. ^1 the "Fermi surface" in LSCO 
evolves from the "Fermi arc" in lightly-doped samples, 
to a hole-like Fermi surface in underdoped samples, and 
to an electron-like Fermi surface in overdoped samples 
(x>0.15). 

The evolution of electronic structure with doping in 
(Ca2-xNa2:)Cu02Cl2 ex hibit s marked resemblance to 
that in (La9-^Sr^)CuOA l)l63l) . As summarized in Fig. 
1451 at low doping, the quasiparticle weight is again con- 
fined to the nodal region and the weight in the quasipar- 
ticle peak, Zgp, increases with increasing doping, consis- 
tent with LSCO. In the Na— CCOC system, recent scan- 
ning tunneling microscopy (STM) work has revealed a 
real spa ce pa ttern of 4ao x 4ao two-dimensional charge 
ordering('l82'). In momentum space, as seen from Fig. 
1461 strong Fermi surface nesting exists in Na— CCOC 
with a nesting vector insensitive to doping close to 2*7r/4 
that may account for the broad near-Ej? spectra in the 
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FIG. 43 Dichotomy between nodal and antinodal excitations 
in LSCO l|l34|) . (a). Spectral weight near a small energy 
window of Fermi level as a function of fc^ and ky for LSCO 
X — 0.063 sample measured at ~20K. The original data was 
taken in the second Brillouin zone and converted into the 
first Brillouin zone and symmetrized under four-fold symme- 
try, (b). Experimental Fermi surface for LSCO x — 0.063 
sample. The black open circles are obtained from the MDC 
peak position at Ep- The solid lines are guides to the eye for 
the measured Fermi surface. The red lines represent the por- 
tion of Fermi surface where one can see quasiparticle peaks. 
The dotted black line represents the antiferromagnetic Bril- 
louin zone boundary; its intersection with the Fermi surface 
gives eight hot spots (solid yellow circles) from {iT,n) magnetic 
excitations. The double-arrow-ended green line represents a 
nesting vector between the antinodal part of the Fermi sur- 
face (c). EDCs on Fermi surface for LSCO x = 0.063 (cl), 
0.09 (c2), and 0.22 (c3) samples. All samples are measured at 
~20 K. The corresponding momentum position is marked in 
the upper inset of each panel. Also included are the spectra 
at (7r,0) points, colored as blue. 



anti-nodal region. In LSCO, neutron scattering has also 
indicated the existence of dynamic stripes ! 183i) . These 
similarities suggest an intrinsic commonality between the 
low-lying excitations across different cuprate families and 
may imply a generic microscopic origin for these essen- 
tial nodal states irrespective of other ordering tenden- 
cies. At very low doping levels, the nodal excitations 
should entirely dominate the transport properties, con- 
sistent with the high-temperature me tallic tendencies ob- 
served in very lightly doped cuprates l|184|) . Thus any mi- 
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FIG. 44 Doping evolution of underlying "Fermi surface" in (La^-TSr^lCuOii ilS' 
a temperature of ~20K. 
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croscopic models of charge ordering must simultaneously 
explain and incorporate the existence of coherent nodal 
states and broad antinodal excitations. 

The nodal-antinodal dichotomy of quasiparticl e dy - 
namics in the normal state also exists in Bi2212 l|l85j) . 
A number of possible mechanisms have been proposed 
to account for the antinodal spectral broadening in the 
normal state. A prime candidate is the (tt^tt) mag netic 
excitations observed in various cuDrates l|79t ISOl: l8ll) . As 
schematically shown in Fig. 143b . this excitation will 
give rise to "hot spots" on the Fermi surface that can 
be connected by (tTjTt) momentum transfer. Electrons 
around these hot spots experience additional scattering 
from the (tTjTt) magnetic scattering. The same mech- 
anism has also been proposed for (Nd2-a:Cea;)Cu04 for 
which the spectral broadening is localized to the expected 
"hot spot" (186). However, in LSCO, the same magnetic 
response, magnetic resonance mode, is not observed. In- 
stead, incommensurate magnetic peaks are observed at 
low energy (below 15 meV)(18J|), which broaden rapidly 
with increasing energy althoug h the magnetic fluctua- 
tion can persist up to 280 meV ()l87l) . Intrigued by the 
fact that the extra broadening sets in when the Fermi 
surface turns from the (7r,0)-(0,7r) diagonal direction to 
the (0,0)-(7r,0) or the (0,0)-(0,7r) direction [Fig. EHb and 
c), an alternative mechanism was proposed(|l3^ in which 
the scattering in question causes a pair of electrons on two 
parallel antinodal segments to be scattered to the oppo- 
site ones (Fig. 143b). In the normal state, this scattering 
can cause a quasiparticle to decay into two quasiparticles 
and one quasihole. The antinodal spectral broadening oc- 
curs as a result of the frequent occurrence of such a decay 
which renders the normal state quasiparticle ill defined. 



Another potential explanation for the broad antinodal 
features may come from models based on the polaron 
picture discussed before l) 11 St Il67l : Il75|) . In such a sce- 
nario, the strong coupling of the electrons to any bosonic 
excitations would result in Z <Cl, and spectral weight 
is transferred to incoherent, multiboson excitations. An 
effective anisotropic coupling could lead to a larger Z 
(weaker coupling) along the nodal direction and a much 
smaller, yet still finite Z, at the antinodes (strong cou- 
pling). In this picture, the antinodal polaron effect in 
LSCO (Fi g. |B ;) |T3|) is much weaker than Na-CCOC 
(Fig. 145b ) flGSJ) if one compares the spectral weight near 
Fermi level around the antinodal region. Regardless of 
the microscopic explanation, the broad and nested antin- 
odal FS segments observed by ARPES are consistent 
with the propensity for 2-dimensional charge ordering in 
the lightly dope d cuprates seen in STM exp eriments on 
Na-CC0C(|l8l and Bi2212(|l8i IHS [190). Further- 
more, an explanation based on an anisotropic coupling 
(coming from either polaron physics or the magnetic res- 
onance) may not be sufficient to cause the 2-dimensional 
charge order; it may be a combination of strong coupling 
and Fermi surface nesting which ultimately stabilizes the 
antinodal charge-ordered state. 



E. Electron-Phonon Coupling and High Temperature 
Superconductivity 

Much of the physics discussed in this review has 
attributed essential features of the ARPES data to 
electron-phonon coupling, and if not to electron-phonon 
coupling alone, to electron-phonon coupling in an an- 
tiferromagnetic background. The question remains as 
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FIG. 45 Dich otom y between nodal and antinodal excitations 
in Na-CCOC ljl63l) . (a). Schematic of the low-lying spectral 
intensity for (Ca2_xNa^)Cu02Cl2 (x=0.10). The hatched 
regions show the nested portions of Fermi surface, and the 
Fermi surface angle is defined in the lower right quadrant, 
(b). EDCs taken at equal increments along the FS contour 
from the nodal direction (top) to the antinodal region (bot- 
tom) for x=0.05, 0.10, and 0.12 at a temperature of 15 K. 
(c). The doping evolution of the low- lying spectral weight 
(circles), along with corresponding data from La2Ca;Sri:Cu04 
(squares), with the error bars representing the uncertainty in 
integrated weight as well as sample-to-sample variations. 




FIG. 46 Doping evolution of "Fermi surface" in Na- 
CC0C( 16iD . (A to C) The momentum distribution of spec- 
tral weight within a ±10-meV window around E_f for x=0.05, 
0.10, and 0.12 in one quadrant of the first Brillouin zone. Data 
were taken at 15 K and symmetrized along the (0,0)- (7r,7r) 
line. The data acquisition range is shown within the black 
lines. The FS contours shown in (D to F) were compiled 
from more than four samples for each composition with dif- 
ferent photon energies and photon polarizations. Data from 
these samples constitute the individual points; the best fit is 
shown as a solid line. The region in which a low-energy peak 
was typically observed is marked by gold circles. The gray 
shaded areas in (E) represent the momentum distribution of 
intensity at Ef±10 meV along the (0,0)-(7r,7r) and (7r,0)-(7r,7r) 
high-symmetry directions. 



to how this electron-phonon coupling can account for 
high-temperature superconductivity with d-wave pairing 
seen in the cuprates. It is often assumed that el-ph cou- 
pling leads to s-wave pairing, and that therefore such a 
mechanism contradicts with the d-wave symmetry of the 
Cooper pairing in the cuprates. Instead, electronic cor- 
relations have been thought to be consistent with d-wave 
pairing. However, while strong electronic correlations 
will suppress the Cooper pair amplitude on the same or- 
bital, and hence induce a d-wave like symmetry, they do 
not tell us much about the explicit pairing mechanism. 
One of the early studies on pos sible phononic mechanisms 
of high Tc suDerconductivitv l)l54|) pointed out that the 
out-of-plane displacement of the oxygen, i.e., the buck- 
ling mode, combined with antifcrromagnetic correlations, 
leads to d^-^-y-^ pairing. Bulut and Scalapino (,191.) stud- 
ied the various phonon modes from the viewpoint of the 
possible pairing force. They found that the interaction 
which becomes more positive as the momentum trans- 
fer increases helps dx2-y2 pairing (the case for buckling 
mode, but not the case for the apical oxygen mode or the 
in-plane breathing mode). 

One can understand the nature of the q momentum de- 
pendence by considering how the phonon couples to the 
electron density. For deformation phonons, the coupling 
is dipolar driven and thus small for small q, the case for 
the breathing modes. This also includes infrared active 



phonons. Yet for Raman active modes, which couple via 
the creation of isotropic and quadrupolar moments, the 
coupling is generally strongest for small q. Specifically 
for the cuprates, such strong fc, q- dependencies occur 
explicitly for c-axis phonons, which include the Raman 
active in-phase buckling Aig, out-of-phasc buckling Big 
and modes involving the apical oxygen Aig. The k mo- 
mentum dependence comes from the phonon eigenvec- 
tors as well as the direction of charge-transfer induced 
by the phonon. For example, for the Big phonon the 
eigenvectors enforce a change of sign when k^ and ky are 
interchanged, a factor ~ cos{kxa) — cos{kya), while for 
the apical charge transfer coupling between Cu and the 
three oxygen orbitals, a factor ^ [cos{kxa) — cos(A:j,a)]^ 
emerges. 



As discussed by Bulut and Scalapino 1|1 9 1[1 among oth- 
ers, the q dependence of phonons can be important 
to give dj.2_j,2 pairing. In particular, if the attractive 
electron-phonon interaction falls off for momentum trans- 
fers q along the diagonal, then conceptually the interac- 
tion is of the same structure as the magnetic pairing from 
antifcrromagnetic spin fluctuations. This type of struc- 
ture occurs for both Big and Aig c-axis Raman-active 
phonons, and thus they contribute to the pairing inter- 
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action in the d-wave channel, parameterized by A^ 

, _ 2j:,^^,dkdk'\gik,k-k')\^6ie,)6iek') 

with the d-wave form factor d^ = [cos(fc2,a) — cos(fcj^a)]/2. 
However, the Aig phonons predominantly contribute to 
the s-pairing channel (replacing dk by 1 in the above 
equation) in the absence of any Coulomb interaction, 
leaving the Big phonon as the largest contributor to d- 
wave pairing, as found in LDA studies ((122) ■ 

However Coulomb interactions change this picture. 
They cannot be neglected since they are necessary to 
screen the long-wavelength nature of isotropic charge 
fluctuations. The screened electron-phonon interaction 
g is of the form 



g{k,q,n) = g{k,q) + 



V{q)ng^i{q,n) 

i-y(q)ni,i(g,r!)' 



(4.13) 



where V{q) — 47re^/g^ is the 3D Coulomb interaction 
and Ha,fc(g, ri) is the frequency-dependent polarizability 
calculated with vertices a, b respectively. Note if g were 
independent of momentum, then the effective electron- 
phonon coupling would be screened by the dielectric func- 
tion e(g, 17) = l—V{q)Ili,i{q,i^). Particularly in the limit 
g — > we recover complete screening and 5 = for 17 = 0, 
restating particle number conservation, while for 12 = 
^ph the renormalized coupling is of order U,ph/^pi- How- 
ever, any fermion fc-dependence of the electron-phonon 
coupling survives screening even at q = as shown by 
Abrikosov and Genkin l,193.) . and the effective charge ver- 
tex in this limit is g{k, g ^ 0) = g{k, g — > 0) — Sg, with 
Sg = {g{k,q -^ 0)), and (• • • ) denotes an average over 
the Fermi surface, defined as 



{A) 



Y.kA{k)5{e{k)) 



(4.14) 



Thus screening removes the constant part of the electron- 
phonon interaction and can highlight the d-wave chan- 
nel. This is important if the bare coupling is highly 
anisotropic with the Fermionic momentum fc, the case 
of the apical oxygen coupling. 

Moreover, the issue of strong local correlations on 
electron-phonon interactions has been recently read- 
dressed by the Hubbard X operator metho d l|90i Il94j) 
and quantum Monte Carlo simulations l|l95j) . Assuming 
no specific phonon and that phonons couple to the on- 
site charge density, i.e., diagonal coupling, these works 
found enhanced forward scattering (i.e., small momen- 
tum transfer), while large momentum transfer process 
were suppressed. Therefore, d^2_y2 pairing can occur by 
el-ph coupling. Furthermore, the vertex correction ex- 
plains the absence of phonon features in the resistivity, 
since the transport relaxation rate contains the factor 
1 — cosO {9: the angle between the initial and final state 
momenta) which reduces the contribution for forward 
scatterings. There has been controversy as to whether 



the vertex correction for the off-diagonal el-ph coupling, 
which modulates the bond, also enhances forward scat- 
tering and suppresses large momentum transfers l|llCl|) . 
The in-plane half-breathing mode, which modulates the 
bond, exh ibits a sharp softening with doping in neutron 
scattering ()l07j) and has been studied in particular. The 
Zhang-Rice singlet couples to the half-breathing mode 
much stronger than estimated in LDA calculations, and 
the vertex correction leads t o an effective attractive in- 
teraction for d^2_y2 pairing l|l96() . On the other hand, 
later analysis (197) has shown that the cancelation of 
terms reduces the off-diagonal coupling, and the diago- 
nal coupling dominates even after the vertex correction 
has been taken into account. In understanding the ef- 
fects of this vertex correction on experimental spectra, 
one should note that the correction works differentl y for 
phononic and electronic self-energies. The sum rules (|198|) 
conclude that the phononic self-energy is reduced by an 
additional factor of x (hole concentration) as compared 
to the electronic self-energy. Intuitively, the difference 
between phononic and electronic self-energies arises be- 
cause a small number of holes cannot influence phonons 
as much as phonons, in which atoms vibrate at every site, 
can influence a single hole. 

In summary, local coulomb repulsion suppresses 
charge density modulations, which in turn decreases the 
strength of the electron-phonon interaction at large mo- 
mentum transfers. This has two effects: first, as a con- 
sequence the contribution of all phonons to the resistiv- 
ity will be reduced by the correlation effect. Second, 
and more relevant to pairing, small q phonons will have 
an accentuating A for d-wave pairing since the coupling 
will decrease faster for large q than without correlations. 
Thus it appears that Coulomb interactions in general 
can have a dramatic impact on electron-phonon driven 
dr^2_y2 pairing. However theoretical developments are 
still needed in order to treat the simultaneous impor- 
tance of strong correlations and electron-phonon cou- 
pling. This is a promising direction for future research. 



V. SUMMARY 

ARPES experiments have been instrumental in iden- 
tifying the electronic structure, observing and detail- 
ing the electron-phonon mode coupling behavior, and 
mapping the doping evolution of the high-Tc cuprates. 
The spectra evolve from the strongly coupled, pola- 
ronic spectra seen in underdoped cuprates to the Migdal- 
Eliashberg like spectra seen in the optimally and over- 
doped cuprates. In addition to the marked doping de- 
pendence, the cuprates exhibit pronounced anisotropy 
with direction in the Brillouin zone: sharp quasiparti- 
cles along the nodal direction that broaden significantly 
in the anti-nodal region of the underdoped cuprates, an 
anisotropic electron-phonon coupling vertex for particu- 
lar modes identified in the optimal and overdoped com- 
pounds, and preferential scattering across the two paral- 
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lei pieces of Fermi surface in the antinodal region for all 
doping levels. This also contributes to the pseudogap ef- 
fect. To the extent that the Migdal-Eliashberg picture 
applies, the spectra of the cuprates bear resemblance 
to that seen in established strongly coupled electron- 
phonon superconductors such as Pb. On the other hand, 
the cuprates deviate from this conventional picture. In 
the underdoped regime, the carriers are best understood 
as small polarons in an antiferromagnetic, highly elec- 
tron correlated background, while the doped compounds 
require an anisotropic electron-phonon vertex to detail 
the prominent mode coupling signatures in the super- 
conducting state. Electronic vertex corrections to the 
electron-phonon coupling furthermore may enhance, and 
for certain phonons, determine, the anisotropy of the 
electron-phonon coupling. A consistent picture emerges 
of the cuprates, combining strong, anisotropic electron- 
phonon coupling, particular phonon modes that could 
give rise to such a coupling, and an electron-electron 
interaction modifying the el-ph vertex. Such a combi- 
nation, albeit with further experimental and theoretical 
effort, may indeed lead to an understanding of the high- 
critical transition temperature with d-wave pairing in the 
cuprate superconductors. 
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